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Description 

Field of the Invention 

5 This invention relates to method and apparatus for reducing or eliminating spurious artifacts tending to 
appear in two-dimensional displays of logging information obtained by moving an array of sensors along a sur- 
face of an earth formation, for example along the wall of a bore-hole made in the course of exploring for pet- 
rochemical substances. 

10 Background of the Invention 

In the exploration of earth formations, it is known to move an array of sensors along an exposed surface 
of the earth formation, typically along the exposed wall of a bore-hole in the earth, to produce data signals 
representative of variations in certain meaningful characteristics of the adjacent earth formation. For example, 

15 this is commonly done in the so-called logging of bore holes, in search of natural deposits of oil. Typically also, 
the data signals produced by the sensors are in electrical form, or are converted thereto. 

A common form of sensor system used for this purpose injects current from an electrode array into the 
adjacent earth formation and measures variations in that current and/or in a voltage related to the current. It 
is also possible to use an array of electrodes to sense variations in a potential known as the spontaneous po- 

20 tential or SP, which appears naturally along the explored wall. In either case, the electrical signals vary in ac- 
cordance with certain characteristics of the adjacent earth formation, and can therefore be used to characterize 
that formation. 

While it is possible to obtain some useful information in such systems which use only a single sensor, sig- 
nificant advantages result from using systems containing an array of multiple sensors to obtain the data signals, 

25 and using the data signals so derived to produce a two-dimensional image corresponding to an area of the 
bore-hole wall extending both along, and at right angles to, the direction of motion of the sensor array. One 
such system, commonly known as a Formation Microscanneror FMS™, conveys to an experienced observer 
a clearer mental picture of the distribution of the measured characteristic over a segment of the bore-hole wall 
than does a single-line graph of data produced by a single sensor. 

30 One way in which to collect data to produce a two-dimensional image of a bore-hole wall would be to use 
a linear array of sensors, i.e., on in which the sensors are all distributed along a straight line, at right angle to 
the direction of motion of the array. In this way, the data collected would span two dimensions: at any given 
time the array of sensors would be collecting data along the azimuth direction (i.e., the direction transverse 
to the array motion along the bore-hole wall) and, over time, the array would be moved to successive different 

35 positions along the bore-hole wall. 

In any sensor system, the individual sensors need to be sufficiently separated one from another to reduce 
the possibility that the sensors will interfere with one another in their operation; for example, sensor electrodes 
will electrically short-circuit to each other if they are placed too close together in the array. However, a dele- 
terious result of such sensor separation in a linear array is the decrease in the resolution of the final image in 

40 the azimuth direction, or an increase in undesirable aliasing problems; in this connection see U.S. Patent No. 
4,567,759 of Ekstrom et al, issued February 4, 1986. 

One way to accommodate the sensor separation requirement and still achieve satisfactory azimuth reso- 
lution is to use a system in which the array of sensors is distributed in a two-dimensional pattern, so that the 
sensors span the azimuth direction with increased resolution (i.e., there are a larger number of sensors per 

45 u nit of distance along the azimuth direction) while still being sufficiently separated one from another. One such 
sensor pattern positions the sensors across the array in a zig-zag pattern, whereby the sensors are very close- 
ly spaced, or even overlapping, along the azimuth direction, but are separated into rows such that adjacent 
azimuth sensors are in different rows. As such a sensor array is pulled along the bore-hole wall, the sensor 
array sweeps out a continuous, or more clearly continuous, area of the bore-hole wall than if a single row of 

50 sensors was used, and yet provides the necessary separation between electrodes. 

At any given time, the array of sensors will be producing data corresponding to different positions on the 
bore-hole wall, both in the azimuth direction and in the logging direction, corresponding to the different posi- 
tions of the sensors distributed on the face of the probe. For simplicity, position along the logging direction 
will be referred to herein as "depth," even though the logging direction is not necessarily limited to being purely 

55 vertical. In order to achieve maximal resolution, these different positions need to be taken into account when 
the data is presented in image form. More particularly, the depths associated with the sensor data from indi- 
vidual sensors need to be adjusted to ensure that all the data are presented for depths corresponding to those 
at which the sensors were when the data was produced. A system of this type using depth adjustment of the 
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data is described, for example, in the above cited Ekstrom et al patent. 

While two-dimensional display systems have proved to be effective for their intended purposes, as with 
nearly all measuring systems their maximum useful sensitivity and accuracy is limited by the presence in the 
data signals of variations which do not represent the characteristics to be detected. Such undesired, false, or 
5 spurious signal components will be designated herein as "noise," to distinguish them from "true" data signals, 
without thereby suggesting that they are necessarily random. 

When the desired, true data signals are large compared to the level of noise, noise is not as much of a 
problem as when the desired, true data signals are relatively small. Forexample, in usual bore-hole exploration, 
a water-based "mud" is used in the bore-hole which helps to support the bore-hole walls physically, and also 
10 provides a medium of rather high electrical conductivity for current flowing between sensors and bore-hole 
wall. With such a high-conductivity mud, it is relatively easy to obtain data signals sufficiently large compared 
to the level of noise to be practical for commercial FMS™ bore-hole exploration. However, it is sometimes de- 
sirable for certain purposes to use an oil- based mud, which has a much lower conductivity than usual water- 
based muds, and this lower conductivity reduces the strength of the desired, true data signals so that they 
15 are more easily obscured by noise signals. Accordingly, in this situation, it is especially desirable to discriminate 
against noise signals and to reduce their effects on the two-dimensional signal display. 

The problem of noise affecting formation logging signals is not unique to the situation described above 
and various approaches have been suggested to deal with the problem. In US 3,521,154, enhancement of log- 
ging signals is attempted by making multiple measurements of the same formation characteristic, in this case 
20 formation resistivity, at a given location in the borehole and combining these signals to enhance the signal rep- 
resenting the formation characteristic. In US 3,457,500, a method is disclosed in which the error (noise) af- 
fecting logging signals is estimated by comparing a sequence of signals in which the error is alternately additive 
and subtractive of the logging signal. The estimated error can then be removed from the logging signal. 

One type of noise which can adversely affect the results of borehole exploration is designated herein as 
25 "one-dimensional" noise, since it varies primarily as a function of the position of the entire sensor array along 
its direction of motion, but is substantially the same for all sensor positions in the array. When such one-di- 
mensional noise occurs in a system in which the sensors are distributed in a two-dimensional array, all of the 
sensors will be affected by the noise at the same time. After the data is depth-adjusted to take into account 
the vertical distribution of the sensors on the probe face, the noisy data will appear as anomalous artifacts in 
30 the final image in the shape of "footprints" of the pattern of sensors in the array. For example, where the sen- 
sors extend across the array in a zigzag pattern, the artifacts are in the form of a corresponding set of zig- 
zags. 

Accordingly, it is an object of the present invention to provide a new and useful signal filtering method and 
apparatus for reducing or eliminating certain artifacts which may appear in an electrically-derived two-dimen- 
35 sional image of an earth formation due to one-dimensional noise signals present in a two-dimensional array 
of logging sensors. 

Another object is to provide a method and apparatus for reducing or eliminating artifacts of zig-zag form 
which tend to be produced in such a two-dimensional image by one-dimensional noise signals from a two- 
dimensional array of logging sensors arranged along zig-zag lines. A further object is to provide a method and 
40 apparatus for producing more accurate, and more readily interpreted, two-dimensional images of certain char- 
acteristics of a borehole wall. 

Summary of the Invention 

45 In accordance with a first aspect of the present invention, there is provided a method for producing a two- 
dimensional image of a characteristic of an area of a wall of an earth formation using an array of sensors com- 
prising a plurality of sets of one or more sensors, some of said sets being at different positions in said array 
as measured along said first direction, and at least one of the sensors of at least one of said sets being dif- 
ferently positioned from all of the sensors in at least one of the other of said sets, as measured along a second 

so direction transverse to said first direction, the method comprising moving the array along said wall in a first 
direction to produce a signal level at each sensor representative of the value of said characteristic in the portion 
of said wall adjacent thereto, sampling said signal levels at all of said sensors of said sets at each of a suc- 
cession of sampling times, to produce successive original frames of sensor signal samples, processing said 
samples and applying said samples to a two-dimensional display device in sets such that said processed sam- 

55 pies of each of said original frames control points in said display which are in a geometric pattern corresponding 
to the geometric pattern of said sets of sensors in said array, characterized by, prior to applying said processed 
samples to said image-display, generating from said samples a correcting signal n, and subtracting it from 
each of said samples, to reduce any artifact in said image produced by noise affecting said sensors, said noise 
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affecting each of said sensors substantially equally during any of said sampling times and producing an artifact 
in said image having a form corresponding to said pattern of said array of sensors, wherein the values of n 
are generated such as to minimize abrupt discontinuities between values of those samples which are produced 
in different original frames by sensors at the same position along said first direction, and wherein the values 
5 of n so subtracted are equal for all samples produced at the same sampling time. 

In accordance with a second aspect of the present invention, there is provided a system for producing a 
two-dimensional image of at least a segment of a surface of an earth-formation wall, comprising (a) a two- 
dimensional array of sensor electrodes and means for moving said array along said wall segment in a first di- 
rection, the center of said electrodes being spaced from each other both along said first direction and along 
10 a second direction at right angles to said first direction, at least some of said electrodes at different positions 
along said first direction being non-aligned with each other along said first direction, (b) sampling means for 
repetitively sampling the signal levels at all of said electrodes to form corresponding original frames of signal 
samples, (c) depth-adjustment means for processing said original frames of signal samples to produce depth- 
adjusted frames of said signal samples, all of said samples in each of said depth-adjusted frames being pro- 
fs duced at the same position along said wall, (d) a two-dimensional display device responsive to said depth- 
adjusted frames of signal samples to produce successive lines in said image in response to corresponding dif- 
ferent ones of said depth-adjusted frames of signal samples, thereby to form a two-dimensional image of said 
wall segment, characterised by (e) filter means for filtering said depth-adjusted signal samples to reduce ar- 
tifacts introduced into said image by noise which affects all of said electrode signal levels substantially equally, 
20 said artifacts corresponding to the patterns of electrodes in the array, said filter means comprising noise signal 
estimating means responsive to said depth-adjusted samplesto generate correcting signals having discrete 
values n(k) identified with each of said depth-adjusted samples, and means for subtracting each value n(k) of 
said correcting signals from the corresponding sample to reduce said artifacts in said image, and wherein said 
correcting signal values are selected to minimize differences between the values of said depth-adjusted signal 
25 samples produced by electrodes adjacent to each other along said second direction, and to distribute said cor- 
recting signals equally oversets of depth-adjusted samples produced by the same original frames of said signal 
samples. 

The objects of the invention are obtained by the provision of a filtering method and apparatus which filter 
the combined noise and true data signals produced by an array of staggered sensors, to reduce or eliminate 
30 undesired artifacts otherwise produced in the two-dimensional image after depth adjustment of the noise and 
data signals from the sensors. This filtering method and apparatus makes use of the facts that the so-called 
one-dimensional noise components from the sensors are generated all at the same time, while data signals 
for any given geological feature at a given depth along the wall of the earth formation are captured at different, 
but known, times by sensors in different rows of the array; that measured data containing both noise and true 
35 data are depth-shifted by the same known amounts prior to application to the image-display device, thereby 
giving the undesired noise artifacts their characteristic shape; andthat geologic features tend to change only 
slightly and/or gradually along the azimuth direction, e.g., in horizontal planes. In performing its function, the 
filter method and apparatus of the invention take advantage of these facts in generating correcting signals n 
corresponding to the noise components which produce the above-described artifacts, these correcting signals 
40 n being subtracted from the measured signals M to arrive at the desired, true data signals S = M - n. The values 
of n so subtracted are equal for all samples produced at the same sampling time. 
More particularly, S(i,j) = M(i j) - n(k), where: 
i is an integer index identifying the depth at which an electrode sample is taken, preferably in units cor- 
responding to the distance the array moves between samplings of the signal levels at the sensors of the array; 
45 j is an integer index identifying the azimuth position of the sensors in the array (e.g., at right angles to 

the direction of motion of the array), numbered from one side to other of the array; 

k is an integer index identifying the depth of the array corresponding to the time when the data was 
collected at position (i,j); S(i,j) represents the true data values in the signal for position (i,j); 
M(ij) represents the measured values of true signal plus noise for position (i,j); and 
50 n(k) is the electronically generated, estimated value for the one-dimensional noise in the frame of data 

collected when the array was at depth k. 

k is preferably represented as i-g(j), where g(j) is a function identifying the relative depth of each sensor 
at azimuth position j, measured with respect to the sensors of the top row, this relative depth preferably being 
an integer multiple of the distance the array moves between samplings, that is, g(j) is preferably in units of i. 
55 As such, k, the depth of the array, is preferably also the depth of the sensors of the top row of the array, and 
g(j) is preferably 0 for those azimuth positions j which correspond to sensors in the top row. 

For example, in a case described hereinafter using four rows of electrodes equally spaced along the di- 
rection of logging motion, and with the array moved by the full inter-row spacing between successive sampling, 
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g(j) may be expressed as R[j], where R[j] is the remainder of the quotient j/4 and j is numbered consecutively 
in the azimuth direction, from one side of the array to the other, starting with an electrode at one end which 
is in the second row, as will appear later herein; in this case k = i-R[j], each value of k identifying a different 
zig-zag set of simultaneously produced data samples among the depth-adjusted samples M(i,j). 
5 Each true signal value S(i,j) is obtained by subtracting from the measured value M(i,j) an estimated noise 

value n(k) = n[i-g(j)] evaluated for the same values of i and j. Since geologic features change only slowly along 
the azimuth direction, the true data signals produced by azimuthally-adjacent sensors should be roughly the 
same, i.e.: 

S(ij) - S(ij + 1)*0. 

10 Using this fact signals having the values n(k) are generated from the M(i,j) signal samples and from the array 
configuration function g(j) in a manner which minimizes the discontinuities, or differences, between the values 
of depth-adjusted samples produced by azimuthally-adjacent sensors. 

This preferably done by generating in a computer the values of n(k) which minimize the image discontinuity 
measure I, where: 

15 

1=11 {S(i,j) - S(i,j+1)) 2 , 

i j 

20 or, substituting M-n for S, 



1 = 21 {K(i,j) - n[i-g<j)] - M(i,j+l> + n(i-<j(j+l)l) 2 . 



The image discontinuity I is minimized when its partial derivatives with respect to the estimated noise n(k) 
are zero. Accordingly, n(k) preferably has values which satisfy the set of linear simultaneous equations ob- 
tained by setting dl/dn(k) = 0, and set forth in detail hereinafter. 
30 The solution n(k) to this set of linear simultaneous equations is expressible through mathematical manip- 
ulation as: 

D 

n<K) = I f'fJc.x) • y(x) , 



where T _1 is the inverse of the matrix T which represents in matrix form the coefficients of the noise terms in 
the above-mentioned set of linear simultaneous equations, x represents the summation parameter and runs 
from 1 to D, where D is the total number of frames in the image, k and x are also the indices of the inverse 

40 matrix T~ 1 , and y is a function of x containing all terms in M resulting from the above-described partial differ- 
entiation of the discontinuity measure I. 

In the preferred embodiment, the computer does not perform all of the many calculations implicit in the 
last equation for n(k), since they are so many and so burdensome. Instead, it is preferred to evaluate n(k) using 
a mathematical sliding window of width w. The smaller the value of w selected, the simpler is the computation 

45 of n(k), although the accuracy of the calculation may decrease when smaller values of w are used. Areasonable 
compromise in many cases is to select w = 11. The preferred form for n(k) is then: 

v v 

n(k) = i/w - z ?- l (k t ,k t ) . y (*♦*.-**> # 

50 k 

where kj and k 2 are the indices of the inverse matrix T-\ and y is now a function of k+k r k 2 containing all terms 
in M resulting from the above-described partial differentiation of the discontinuity measure I. 

55 

Brief Description of Figures 

These and other objects and features of the invention will be more readily understood from a consideration 
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of the following detailed description, taken with the accompanying drawings, in which: 

Figure 1 is a schematic vertical section showing a bore hole and apparatus for conducting logging of the 
bore-hole wall, in which the present invention is useful; 

Figure 1 A is an enlarged fragmentary view showing the part of the system of Figure 1 which carries the 
5 logging array; 

Figure 1B is a fragmentary perspective view showing a probe pad which may be used in the apparatus of 
Figure 1A; 

Figure 2 is a block diagram of the electrical elements of the system of Figure 1, including apparatus op- 
erating in accordance with the invention; 
10 Figures 3, 4, and 5 are tables and diagrams to which references will be made in explaining the nature of 
the problem which the invention solves; 

Figure 6 is a schematic diagram of a sensor array, showing a representative electrode configuration and 
specific array dimensions usable in connection with one application of the invention. 
Figure 7 is a table of hypothetical, assumed values of electrode signal samples at various depths of the 
15 array, to which reference will be made in explaining the manner in which the invention in one of its forms 
may be applied; 

Figure 8 is a flow diagram illustrating the basic steps of a preferred procedure for practicing of the invention 
in one of its forms; Figure 9 is a flow diagram illustrating in somewhat more detail one manner of practicing 
the method of the invention in one of its forms; 
20 Figure 10 is a flow diagram illustrating in even more detail a preferred way of practicing the method of the 
invention; 

Figure 11 is a block diagram illustrating in apparatus form a preferred embodiment of the invention; 
Figure 12 shows pairs of logging images obtained using a water-based mud, one pair before and one pair 
after insertion of the filter of the invention into the logging system, respectively. 

25 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Overall System 

30 Referring now to the specific preferred embodiments of the invention shown in the drawings by way of 
example only, Figure 1 illustrates one known type of system to which the invention is applicable. In this system 
alogging tool 10 carries a set of four pads such as 12 which are urged outwardly against the interior wall 14 
of a bore-hole 16 extending vertically into the earth 17. The tool and the pads are pulled upwardly by a pulling 
cable 22, which extends around appropriate sheaves such as 23, 24, supported by tower 25, to a winch 26 on 

35 a vehicle 28 having an engine 29 for operating the winch, preferably to pull the tool and pads upwardly at a 
uniform velocity. The pulling cable preferably is provided with an outer metallic armor, and contains wires ena- 
bling two-way electrical communication between the vehicle and the tool in the bore hole. 

Each pad 12 carries a two-dimensional array of sensors such as 30, which may be simple recessed elec- 
trodes electrically insulated from each other, positioned to be moved along the bore-hole wall as the pads are 

40 pulled upwardly in contact with the walls. In the case of the system known as the Formation Micros ensor 
(FMS™), a voltage is preferably applied between each electrode and a remote reference electrode, for example 
the conductive armor surrounding the pulling cable, to inject a current into the adjacent earth; in a system for 
logging spontaneous potential SP, no such voltage need be applied. 

Referring particularly to Figure 2, in one simple example the pad 12 may contain 4 rows of electrodes such 

45 as 30, the top row containing 6 electrodes while rows 2 through 4 contain 7 electrodes each, for a total of 27 
electrodes. The rows are staggered with respect to each other azimuthally so as to lie along a sawtooth path, 
as suggested by the schematic sawtooth lines joining the electrodes in Figure 2, where "azimuth" designates 
a direction at right angles to the logging direction, i.e., the horizontal direction for the visual vertical pull. As 
the pad is pulled upwardly, the top row of electrodes traverses six azimuthally spaced-apart vertical lines on 

so the bore-hole wall, while each of the other three rows of electrodes traverses seven different vertical lines, 
azimuthally spaced from each other and from the first row of electrodes, so as to "fill in" the spaces between 
the electrodes of the top row and therefore scan the entire wall segment traversed by the array during the pull- 
ing operation. It will be understood that this resolution of scanning cannot be accomplished by using a single 
row of electrodes, since the electrodes would be electrically short-circuited to each other if placed close enough 

55 together to overlap azimuthally. 

The signals at each of the electrodes are supplied to a sampling circuit, in this example represented as 
sample-and-hold circuit in the form of a frame-grabber 36, which is actuated at predetermined intervals to cap- 
ture and store successive frames of data signal values from all electrodes as the array is being pulled upward. 
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Thus, successive frames of 27 electrode signal values are stored, for example at increments of 0.1 inch of 
motion of the array. The data signals captured in each frame by the frame grabber are then sampled succes- 
sively by MUX 38, which typically also time-multiplexes other conventional data produced by a source 40 of 
other down-hole data, such as accelerometer data useful in connection with determining accurately the depth 

5 of the tool at any time. While the use of a frame grabber is preferred from the viewpoint of obtaining data from 
the array of electrodes at exactly the same times in each frame, it involves some degree of complexity in the 
down-hole apparatus; accordingly, when simplicity is economically important it is preferred to sample the elec- 
trodes sequentially with the MUX during each frame and to compensate for the differences in electrode sam- 
pling times in the up- hole apparatus. 

10 The sampled data from MUX 38 are changed to digital form by digitizer and AID converter 41, and sent 
up- hole to the surface vehicle 28 (Figure 1) over a line within the pulling cable 22. In an electronics unit in the 
vehicle, a conventional digital depth-signal generator 44 produces depth data accurately representing the 
depth of the array at which each frame of electrode signal data is produced, and a signal combiner 48 asso- 
ciated each digital electrode signal sample with a depth number i defining the depth at which the sample was 

15 produced, and with an electrode number j identifying the electrode which produced it. These combined raw 
data are preferably stored, as in a tape recorder 50 or in any other suitable memory device, until the data are 
to be used, whether at the vehicle or at any other convenient data-analyzing station. 

In the Figure 2 system, the tape from the tape recorder is later read in tape reader 52, and the resultant 
data fed to a minicomputer such as VAX computer 54 for example, wherein each sample and its associated 

20 identifying values of i and j are stored in a digital memory 56 and thereafter processed in the signal processor 
60, in the minimal discontinuity filter 62 described below, and in any desired conventional additional signal 
processor 66, prior to application to optical printer 68, which produces a paper output displaying a two-dimen- 
sional image of the segment of the bore-hole wall scanned by the electrode array during the logging pull; the 
printer may also provide individual waveforms of the signals produced at the respective electrodes. A switch 

25 69 is also shown to indicate schematically that the minimal discontinuity filter can be inserted into the system 
as described above, or bypassed by the switch when the operation of the filter is not desired. 

As shown, the mini-computer 54 provides within it the signal processor 60, and also the minimal discon- 
tinuity filter 62 of this invention through which the depth-adjusted electrode signals are passed prior to their 
application to the optical printer 68 where the final two-dimensional image is formed. It will be understood that 

30 while the depth adjustment, filtering, and other signal processing are usually performed by programming of 
the computer 54, they can be performed by appropriate hardware. 

Problem Presented 

35 The nature of the problem which thef ilterof this invention attacks will become clearer from a consideration 
of the schematic and idealized Figures 3 through 5. In Figure 3, depths of the electrode array 12 below the 
surface of the earth are indicated in the column at the left by D1 through D16; the next column to the right 
indicates hypothetical "true 11 signal values assumed to be produced at the electrodes, at corresponding depths, 
in the absence of noise. 

40 To the right of the column of hypothetical, noise-free electrode voltages is shown schematically an array 
of 4 rows of electrodes arranged in a sawtooth configuration, in 13 successively higher positions during a log- 
ging pull. Only two electrodes per horizontal row are shown, rather than the six or seven shown in the array 
12 of Figure 2, and the array is shown at sampling positions spaced vertically from each other by the entire 
distance between rows, in the interest of compactness of representation. The eight electrodes are identified 

45 by azimuth parameter j with values I through VIII. 

More particularly, at time A the staggered array of eight electrodes is in the lowermost position shown, 
and each dot marked A represents an electrode of that array which produces a data signal corresponding to 
a characteristic of the earth foundation at the depth of the electrode at that time. The line segments joining 
the electrodes A are schematic, and indicate that the electrodes so joined are ail part of the same array, and 

so are disposed in staggered positions along a sawtooth line, that is, no two electrodes are aligned with each 
other along the direction of logging motion. It is assumed that the frame grabber of Figure 2 is actuated at the 
times A,B,C,D, etc., so that the signal values, at the electrode positions shown, are sampled, digitized, and 
sent up- hole in digital form, one frame for each position of the array. 

The next-higher position of the array above position A, for which the next frame of electrode data is cap- 

55 tured by the frame grabber, is designated by B in Figure 3, and so on for successively higher positions of the 
array designated as C through M. 

From Figure 3 it will be seen that to obtain eight electrode signals which correspond to the same depth, it 
is necessary to select samples produced during different frames, at different times; for example at depth D1 2, 
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the electrodes providing the information are, reading from left to right, those from frame times E, D, C, B, E, 
D, C and B. 

This relation is depicted in detail in Figure 4, wherein the first, or left-most, column lists the frame sample 
times of Figure 3, and the rows of numbers in columns I through VIII show the signal values on each of the 

5 electrodes at the times indicated in the first column corresponding to the hypothetical true values in the second 
column of Figure 3. The frames of raw data samples shown in Figure 4 are sent up-hole in serial bit form, for 
example, first Frame A from left to right, then Frame B from left to right, and so on. Since the data collected 
atthe same time correspond to positions atfour different depths within the bore hole, the data need to be depth- 
adjusted in order to create a proper image of the bore-hole wall. 

10 For example, to create the portion of the image corresponding to depth D1 2, it is necessary to access data 
collected during sampling times B, C, D and E. Specifically, depth D12 is characterized by selecting the fol- 
lowing electrode-sample time combinations from the table of Figure 4: l-E, ll-D, lll-C, IV-B , V-E, Vl-D, Vll-C, 
and Vlll-B, yielding the data presented for depth D12 in the table of Figure 5. Since the relative positions of 
these samples are known from the geometry of the array, the signal processor in the up-hole electronic is easily 

15 programmed to accomplish such appropriate selection of samples to constitute properly each horizontal line 
of the final two-dimensional image. 

The problem produced by so-called one-dimensional noise, i.e., noise which is a function of time but is 
the same for all electrodes at a given time, is illustrated qualitatively for sample time H in Figures 3 and 4, and 
in the portion of Figure 5 marked by the sawtooth-shaped line. Here it is assumed that at time H only, a strong 

20 noise signal occurs, increasing the signal levels on all electrodes by 40 units. More particularly, as shown in 
Figure 3, all eight electrodes at time H, although at different depths, carry the noise of 40 units superposed 
on their true signal values. As shown in Figure 4 at time H, this frame when multiplexed is made up entirely 
of samples of a high, anomalous value which, if applied directly to the image-display apparatus after depth 
adjustment, will produce a strong sawtooth pattern, as depicted by the dashed line in Figure 5. 

25 More particularly, as shown in Figure 5, depth-adjustment serves to group the data samples into the proper 
depth sets as desired, but the data from the noise-augmented frame at time H are disposed along the zig-zag 
path shown in Figure 5, and specifically along a sawtooth locus corresponding to the positions of the electrodes 
in the array, so as to constitute, in effect, the "footprint" of the array of electrodes. It is this spurious sawtooth 
artifact which the filter of this invention removes or greatly reduces in strength, by generating from the meas- 

30 ured signals certain estimated values for the one-dimensional noise which, when subtracted from the depth- 
adjusted samples, in effect cancel out the undesired, spurious, one-dimensional noise artifacts. 

Analysis for the Minimal Discontinuity Filter 

35 The minimal discontinuity filter of this invention involves the following five aspects of the bore-hole logging 
apparatus and procedure for which this filter is designed: 

1) The configuration of the array used to measure the physical features of the bore hole; 

2) The manner in which data is collected using that array; 

3) The effect that the one-dimensional noise has on the data collected using that array; 
40 4) The nature of that noise for which this minimal discontinuity filter is designed; and 

5) The nature of the physical features of the bore hole that are being measured by that array. 

In this discussion, for the purpose of simplicity of analysis, the usual case will be assumed in which the 
bore hole is vertical. Without limiting the scope of the invention, Figure 6 illustrates various of the parameters 
referred to in the following analysis, for the 27-electrode array shown, for the k th position of the array; the broken 
45 lines indicate the position of the array when the next frame of samples is taken at depth k-1. 

1. Configuration of the Array 

The probe used to measure the physical features of the bore hole includes an array of sensors (e.g., elec- 
so trodes) which are oriented in some known, fixed two-dimensional pattern on the face of the probe. The integer 
parameter j is used to identify the azimuth location of each of the sensors of the array, where the I th and (j+l) 1 * 
sensors are next to each other in terms of azimuth (i.e., horizontal) direction. 

The function gQ) defines the vertical location of the sensors relative to some vertical reference position 
on the probe face as a function of the parameter j. For example, if the vertical reference position is the top- 
55 most sensor in the array, then g(j) is the positive ordinal distance the j* sensor is below the sensor (or sensors) 
at the top of the array of sensors on the probe face, where "top" means highest in the vertical direction. 
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2. Data Collection Using the Bore-Hole Logging Probe 



The bore hole is logged by slowly pulling the probe up out of the hole with the sensors of the probe face 
exposed to the hole surface. As the probe is pulled up the hole, data is collected periodically to sample the 
5 physical features of the hole at fixed depth intervals. If the probe is pulled up the hole at a uniform rate, then 
sampling the data at a fixed time interval will result in the sampling of data at constant depth intervals, repre- 
sented by the constant d. 

If the pattern of sensors on the probe face and the depth sampling interval are selected properly, then, at 
any given time, each sensor will be scanning the bore-hole wall at a unique azimuth position and, except for 
10 the leading or top row of sensors, at a depth that was characterized previously by other sensors at other azimuth 
positions during other data sampling times. The distance (called a) between the centers of adjacent rows of 
sensors is preferably an integral multiple, including 1, of the distance d which the probe travels between suc- 
cessive samplings, so as to assure that data will be collected for all of the azimuth positions at each depth 
position. 

is However, as a practical matter it is not possible to move the probe at a perfectly uniform speed, primarily 
because of variable sticking of the probe along the bore-hole wall combined with the resultant elastic stretching 
and re-contracting of the pulling line; thus even if the rate of pulling of the surface end of the pulling line is 
constant, in general the speed of the probe will be subject to variation. Nevertheless, using techniques now 
known and used in the art, later processing in signal processor 60 produces signals representative of the value 

20 of the logging signal at intervals a along the bore-hole wall, as desired. This is accomplished by using the ac- 
celerometer data to obtain accurate information as to the exact position of the array of electrodes at all times, 
and specifically at the times of the samplings of the electrode signal volume. By interpolating between the 
depth and sample-time data, data providing continuous information as to electrode signal values is generated 
for all depths; by resampling the latter data at the times corresponding to when the array is in the successive 

25 positions separated by the distance a, the desired data is obtained. This data interpolating and resampling is 
done in signal processor 60. For further details of such data processing operations, reference is made to the 
above-cited U.S. Patent No. 4,567,759 and to U.S. Patent No. 4,468,623 of Gianzero et al, issued August 28, 
1984. Accordingly, when reference is made herein to sampling the electrode signal values at certain times or 
at certain depth intervals, this should be taken to include the deriving of comparable samples by the signal 

30 processor 60 from raw down-hole data containing samples which were not necessarily taken at equal depth 
intervals. 

Using the integer parameter k to identify the depth within the bore hole of the vertical reference position 
of the sensor array at any given sampling time, and the integer parameter i to identify the depth within the 
hole of any of the particular sensors at that same time, then the following relation is true: 
35 k = i - g(j), (Eq C1) 

which can be rewritten as: 

i = k * g(j), (Eq C2) 

where both i and k are parameters which correspond to integral numbers of sampling intervals d t and both in- 
crease as depth increases. Equations (C1) and (C2) express the fact that the function g(j) defines the vertical 

40 location of the j" 1 sensor of the array relative to the sensor reference position. If the sensor reference position 
coincides with the top row of sensors in the array and if the j"> sensor is in the top row, then g(j) will be 0 for 
that value of j, and i, the depth of that sensor, will be equal to k, the depth of the top row of sensors. 

The parameter k can also be used to identify the depth of the array itself. When the array is at some depth 
k, the depth i of any sensor in the array can be found by substituting that value k and the appropriate j value 

45 for that sensor into Equation (C2). In this sense, at any given time, all of the sensors in the array are associated 
with the same value k. If the i,j values for each of the sensors are substituted into Equation (C1) for a given 
time, the resulting k values will all be the same, and will correspond to the k value representing the depth of 
the array at that time. Assuming that the vertical reference position coincides with the top row of sensors, the 
k"> data frame is the frame of sensor data collected when the top row of sensors is at depth k (even though 

so not all of the sensors in the entire array are at that depth at that time). 

3. Effect of the One-Dimensional Noise on Measured Signals 

By "one-dimensional" noise is meant noise which is a function of time, but which affects all sensors sub- 
55 stantially equally. This one-dimensional noise affects the measure data signals substantially linearly. That is, 
the measured signal M(i,j) (i.e., the signal measured by sensor j at depth i) can be expressed as: 

M(ij) = S(i,j) + N(i,j), (EqC3) 

where 
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S(i,j) = the ideal or "true" signal corresponding to the characteristic of the bore-hole wall, measured at 
azimuth position j and at depth i, and 

N(i,j) = the noise at sensor j when it is at depth i. Solving for the ideal signal S(i,j), Equation (C3) becomes: 

S(ij) = M(ij) - N(ij). (EqC4) 

4. Nature of the One-Dimensional Noise 



The one-dimensional noise, which this minimal discontinuity filter is designed to reduce, when it does oc- 
cur, affects all of the sensors of the array at the same time and to substantially the same extent. This means 
10 that when the noise occurs, all of the sensors will have the same noise value added to them in accordance 
with Equation (C3). 

Assume an array with E sensors whose relative positions on the probe face are defined by the function 
g(j). If noise occurs during data frame k, then the noise signals N(i,j) for all of the array sensors will have the 
same value, call it n(k), the noise associated with data frame k, where k is related to i and j by Equations (C1) 
15 and (C2). Thus: 

N(i.j) = n[i-gO)] = n(k) (Eq C5) 
for all i,j values associated with the E sensors of data frame k. 



5. Nature of the Physical Features of the Bore Hole 



The physical features of the bore hole that are being measured by the probe typically vary very slowly 
along the azimuth direction. This can be represented mathematically as: 

S(iJ) - S(i,j + 1)«0, (Eq C6) 
where S(i,j) and S(i,j+1) are the ideal signals corresponding to two adjacent azimuth locations (j and j+1) at 
25 the same depth i within the bore hole. 

Derivation of the General Form of the Minimal Discontinuity Filter 

Since the difference between values of adjacent pixels within the image is ideally near zero, as expressed 
30 in Equation (C6), the noise N(iJ) associated with a particular location (i,j) (or image pixel) can be estimated 
by minimizing these differences, also called discontinuities, over the entire collection of data, where each dis- 
continuity is of the form of the left side of Equation (C6). One way to express the total discontinuity over the 
entire image is by the sum of the squares of all the discontinuities over the entire image. Mathematically, this 
is expressed as: 



D E-l 

I = Z Z (S(i,j) - S(i,j+l)r, (Eq C7) 

i=X j=l 



where 

I = the measure of the total discontinuity of the entire image, 

D = the number of depth samplings that composes the entire image, 

and 

45 E = the number of sensors in the array used to obtain the data. 

Substituting Equation (C4) into Equation (C7) yields: 



50 I = I V <H<i,j) - M(i,j) - M<i,j+U + N(i,j + 1)> 2 <*J C8> 
i=l j=l 

Using the fact that the noise is the same for each sensor at a given time, the identity of Equation (C5) is 
used to rewrite Equation (C8) as: 

55 

I = z E t l <H<i,j> - n[i-g(j)J - M(i,j+i> + n[i-g<j+i)]>* ... 
i=i j=i (Eq C " 
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Estimates of the noise n[i-g(j)] can be obtained by minimizing the total discontinuity I. This can be accomplished 
mathematically by setting equal to zero the partial derivative of I with respect to the estimated noise n(k) of 
the generic data frame k, and then solving for the estimated noise, viz: 

Ll- = I 2 <M(i,j) - a[i-g<j>) - Md.j+ll ♦ n[i-g<j+l>1> 

<>M<i,j> 3n[i-g<j>1 ^M{i,j+D ^Z^lttil 



. + — - ) = o 

(Eq CIO) 



but: 

15 = m ^i + ^ = 0. for all values of i, j, and k (Eq C11) 

dn(k) dn(k) J 

an an"(k) @l = °* f ° r ' ' 9(i) n0t = k (Eq C12) 

3DDzM as i (for |.fl fl , Bk (Eq C13) 
20 5n(k) 

d n [ '-9Q + 1 )l = o,fori-g(j + 1)not = k (Eq C14) 
dn(k) 

Miz^Dl = i,f 0 ri-g(j + 1) = k (EqC15) 

25 Substituting Equation (C11) into Equation (C10) and dividing both sides of the equation by 2 yields: 



D E-i 

I I (K(i,j) - n(i-g(j)] - K(i,j+1> + n[i-g(j+l)]) 
i=l j=l 



<)n{i-g<j) ] ^B[i-g(j+l) ] 

+ ) - o (Hq CIS) 



As expressed in Equation (C12) to C(1 5), all the terms in the double summation of Equation (C16) will be zero, 
except for those for which: 

i - g(j) = k (Eq C17) 

or 

40 i-g(j + 1) = k. (EqC18) 

The conditions of Equations (C17) and (C18) can be rewritten: 

i = k + g(j) (EqC19) 

and 

i = k + g(j + 1). (EqC20) 
45 Applying the conditions of Equations (C19) and (C20) to the summation of Equation (C16) by first retaining 
only the terms which satisfy Equation (C1 9) by substituting k+g(j) for i, and then retaining only the terms which 
satisfy Equation (C20) by substituting k+g(j+1) for i, yields: 



50 E-l 

0 = Z <-M(X+g(j) ,j) + n[k+9(j)-9(j>] 

+ M(K+g(j) , j+i) - n[*+g(j)-g(j+l>] 



55 



+ M<)c+g(j+i) ,j) - a[X+g<j+l)-g<j)] 

• H<*+g(j+l) + n{*+g<j*l)-g(j+l>]} (Bq C21) 
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Recognizing that: 

n[k + g(j) - g(D] = n[k + g(j + 1) - g(j + 1)] = n(k), (Eq C22) 
Equation (C21) becomes: 



E-l 

0 = Z (-M(X+g(j) ,j) + nOO + K(X+g(j) ,j+l> - n[*+g< j) -g< ) 
j=l 

+ M(k+g<j+i) , j) - n[k+g(j+l)-g<j)] - M()c+g ( j + l) , + n(k)} 

(Eq C23) 

Rearranging Equation (C23) so that all the estimated noise terms n are on one side and all the measured signal 
M terms are on the other yields: 

15 

E-l 

I <2nU) - n(*+g<j)-g<j+x>] - n[*+g< j+l)-g< j ) ] > = y<k>, 
j=l (Eq C24) 



20 



where 



25 E-l 

y(X) = Z <H()c+g<j) , j) - M(X+g(j) ,j+l) - M(*+g< , j ) + j=l 

M(k+g(j+l) , j+lj ) (Eq C25) 

30 Equation (C24) expresses the relationship between the measured data (as represented by the value of 
y(k) defined by Equation (C25)) and the estimated noise (as represented by the left side of Equation (C24)) 
as a function of the data frame parameter k. Since this relationship is true for all the data collected, that is for 
all values of k, then the relationship of Equation (C24) is the characteristic equation for a set of k simultaneous 
linear equations which can be rewritten in matrix form to cover all the data collected: 

35 y = T • n, (Eq C26) 

where: 



40 



45 



X = 



/ yd) 
y(2) 
y(3) 

y<x-u 
y(*> 
y<*+i) 



50 



I 

I 

I 

I 

I 

I 
I 

/ 



(Eq C27} 



55 



12 



EP0 426 564 B1 



10 



/ nil) \ 

i i 

! n(2) ! 

! n(3) ! 

a = | S j (Eq C28) 

j n(lc-l) j 

I nik) ! 



15 

and 

T = the (D x D) matrix which represents in matrix form the coefficients of the estimated noise terms in the 
set of linear equations defined by the relationship of Equation (C24) and D is the total number of data frames 
20 in the entire image. The estimated noise terms are found by expanding the summation of the left hand side 
of Equation (C24) and reducing by combining similar terms. 

Equation (C26) can be inverted to solve for the estimated noise matrix n: 

n = T- 1 • y , (EqC29) 

where 

25 T" 1 = the inverse of matrix T. 

In principle, a value of n(k) can be obtained for each value of kforvalues of k from 1 to D by solving Equation 
(C29) using the inverse of the entire (D x D) matrix T, such that 

n(k) = T- 1 (k,1)y(1) + T- 1 (k,2)y(2) + T- 1 (k,3) y(3) + ... + T~ 1 (k,k- 1) y(k- 1) + T-i(k,k)y(k) + 
T- 1 (k,k + 1) y(k + 1) + ... + T- 1 (k,D) y(D). (Eq C30) 
30 Equation (C30) can be rewritten in summation form as: 



n(k) 



35 



D 

= r 



t • y<x) 



(Eq C30a) 



where x is the summation parameter which runs from 1 to D. 

However, since k may typically vary from 1 to say 10,000, the process of inverting the (10,000 x 10,000) 
element matrix T is computationally impractical. In practice, it is desirable to select a sliding window of finite 
40 length to meet computational requirements. Equation (C29) is then solved repeatedly within the sliding window. 
If the neighboring windows are shifted from each other by a single depth sample, then the estimated noise 
n(k) averaged over all the overlapping windows is given by a linear filter with depth-independent coefficients: 



45 v v 

n<k) = l/w Z Z {?' l i)c u X t ) - Y(k+k.-k,)}, (Bq C31) 

fc T =l k^l 

where 

so k = parameter identifying the frame of measured data, 

n(k) = estimated noise associated with the k" 1 frame of data, 

w = the width of the filter window, 

k 1f k 2 = the summation variables, 

T -1 ( k i. k 2) = the (k 1t k2) element of the matrix T- 1 , 
55 T-i = the inverse of the banded symmetric Toeplitz matrix T, 

T = the (w x w) matrix which represents in matrix form the coefficients of the estimated noise terms in 
the set of linear equations which relate estimated noise to the measure signals, and 

y(k+k r k 2 ) = function of the measured data given in Equation (C25). 
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To demonstrate the implementation of these mathematical steps, the next section presents the derivation 
of the minimal discontinuity filter for a selected electrode array according to Figure 6, used in a particular fash- 
ion, and subject to specified hypothetical signal and noise levels as illustrated in Figure 7. 

5 Derivation of a Particular Minimal Discontinuity Filter (1 inch = 2*54cm) 

One particular example of an electrode array that can be used in the logging of bore holes is shown in 
Figure 6. It contains 27 electrodes arranged in four horizontal rows, designated as Rows 1, 2, 3, and 4, from 
the top downward. The centers of the electrodes are spaced vertically from each other by a distance a = 0.4 

10 inches and horizontally by about 0.1 inch. The 27 electrodes, which are circular and about 0.2 inches (5 milli- 
meters) in diameter, are identified by the azimuth position j, where the first electrode, the (j=1)-electrode, is 
the right-most electrode in the second row (as viewed in Figure 6). The (j=2)-electrode, which lies in the third 
row, is displaced in azimuth direction from the (j=1)-electrode by 0.1 inches. This pattern continues for all 27 
electrodes, such that Row 1 contains 6 electrodes 0=4,8, 12,1 6,20,24), Row 2 contains 7 electrodes 

15 0=1,5,9,1 3,1 7,21 ,25), Row 3 contains 7 electrodes 0=2,6,10,14,18,22,26), and Row 4 contains 7 electrodes 
0=3,7,11,15,19,23,27). 

Data is collected by measuring the signal from all of the 27 electrodes at a given time. Thus, a given data 
frame (identified by the letter k) of measurements from all 27 electrodes contains data for four different depths 
(identified by the parameter i) within the hole. The signals measured with the 6 electrodes of Row 1 will cor- 
20 respond to a particular depth in the hole of, for example, 8.4 inches. The 7 signals from Row 2 will correspond 
to a depth in the hole 0.4 inches (=a) below those recorded by the Row 1 electrodes. Similarly, the Row 3 
electrodes will be measuring data for a depth 0.4 inches below those of Row 2; and Row 4 t 0.4 inches below 
those. 

If the signals are recorded at particular depth intervals d as the probe is pulled up the hole, then the data 

25 can be collected so as to maximize the characterization of each depth selected. Using the array depicted in 
Figure 6, this can be achieved, for example, by recording data every time the probe moves 0.4 inches up the 
hole. For convenience in the present calculation, it is assumed that the depth interval d between successive 
samplings equals the inter-row distance a. 

For example, assume that Row 1 of the array is at hole depth 8.4 inches when the k th data frame is collected. 

30 This depth can be identified by the hole depth parameter i. As explained above, Row 2 will be at depth i+1 
which is 8.4"+ a = 8.4"+0.4" = 8.8", Row 3 at i+2 which is 8.4"+ 2A = 8.4"+0.8" = 9.2", and Row 4 at i+3 which 
is 8.4"+ 3a = 8.4"+ 1.2" = 9.6". When the next data frame (k+1) is collected, after the probe moves d = 0.4 inches 
up the hole, then the Row 1 data will corresponds to depth i-1, which is 8.4"- a = 8.4"-0.4" = 8.0"; Row 2 data 
will correspond to i, which is 8.4"; Row 3 data will correspond to i+1, which is 8.8"; and Row 4 data will corre- 

35 spond to i+2, which is 9.2". One can see that as the probe is pulled up the hole, the selection of 0.4 inches as 
the sampling depth interval d results in eventual measurements at each selected depth by all 27 electrodes, 
albeit at 4 different times, thereby maximizing the azimuthal characterization at each selected depth. 

Best characterization in the logging direction is achieved when the distance between electrode rows a is 
an integer multiple of the sampling depth interval d. Thus, for the array in Figure 6, where the distance between 

40 electrode rows is a = 0.4 inches, best depth characterization will be achieved using sampling depth intervals 
d of 0.4 inches, or 0.2 inches, or 0.1333 inches, or 0.1 inches, or 0.08 inches, and so on. Given a known rate 
of upward pulling of the probe, the time interval between successive data frames can be selected to choose 
the desired sampling depth interval d. However, as pointed out above, by extrapolating from the raw data and 
resampling at the desired depth intervals, the signal processor enables selection of a desired depth interval 

45 regardless of pulling rate. 

The one-dimensional noise which the minimal discontinuity filter of this invention addresses is such that 
all 27 electrodes are affected at the same time by the same amount This means that the noise at electrode 
j=1 is the same as the noise at electrode j=2 and electrode j=3, and so on to electrode j=27 for the frame of 
data collected at that same time. If data frame k is collected when Row 1 is at hole depth i = 8.4 inches, then 

50 electrodes j=4,8,1 2,1 6,20,24 will be at hole depth i, electrodes j=1 ,5,9,1 3,17,21,25 will be at hole depth i+1, 
electrodes j=2,6,1 0,1 4, 18,22,26 at hole depth i+2 t and electrodes j=3, 7, 11, 15, 19,23, 27 at hole depth i+3. Any 
noise N(i+1 ,1) (that is, noise at electrode j=1 at depth i+1) will be identical to the noise N(i+2,2) (that is, noise 
at electrode j=2 which is at depth i+2) when Row 1 is at depth i, and so on for the rest of the electrodes. This 
can be expressed mathematically as: 

55 
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N(i+l,l) = N<i+2,2) = N<i+3,3) = N(i,4) = 

H(i+1,5) = N(i+2,6) = H<i+3,7) = N(i,8) = 

N<i+1,9) = NCi+2,10) = N(i+3,11) = N<i,12) = 

5 H<i+l,13) = NU+2,14) = N(i+3,15) = N(i,16) = (Bq A6) 

N<i+l,17) = N(i+2,18) = N(i+3,19) = N(i,20) = 

NU+1,21) = N<i+2,22) = N(i+3,23) = N(i,24) = 

N(i+1,25) = N<i+2,26) = N(i+3,27) 

10 The fact that all 27 noise signals are the same for this one-dimensional noise can be taken advantage of 
by referring to one noise value n(k), where k identifies the data frame which was subjected to that noise level. 
The assumption for this example of the minimal discontinuity filter is that data will be sampled at intervals of 
d = 0.4 inches, so that the probe moves by one electrode row interval (a) between successive frames of meas- 
ured data. This means that the geometry of the electrodes within the array depicted in Figure 6 can be used 

15 to express the data frame parameter k in terms of the hole depth parameter i and the electrode parameter j, 
such that: 

k = i - g(j) = i - RD1. (Eq A7) 

where 

g(j) = the function which defines the vertical location of the array electrodes relative to the top-most 
20 electrode as a function of the parameter j, and 

R[j] = mod (j,4) = the remainder of j divided by 4. (Eq A8) 

So, 

y(X> = 2n<X) - n<X+R[ 1]-R[ 1+1]) - n(X+R[ 1+1]-R( 1]) 
25 + 2n(k) - n(X+R[ 2]-R[ 2+1]) - n(X+R[ 2+l]-R[ 2]) 

+ 2n<X) - n(X+R[ 3]-R[ 3+1]) - n(X+R[ 3+l]-R( 3]) 
+ 2n(X) - n(X+R[ 4]-R( 4+1]) - n(X+R[ 4+l]-R[ 4]) 

(continued) 

30 Now that it is established that, for the particular array and depth sampling interval of this example, the func- 
tion g(j) is equal to R[j] and the number of array sensors is 27, Equations (C24) and (C25) can be rewritten for 
this example, such that: 



35 



40 



45 



50 



£(2a(X) - n(X+R[j]-R[j+l]) - n(X+R[j+l]-R[j] ) } = y(X),(EqA24) 

where 

26 

y(X) = E {M(X+R[j],j) - M(X+R[j],j+l) - H<X+R[j+l],j) + 
j=l 

M<X+R[j+l] , j+1) ) (Bq A25) 

Expanding Equation (A24) yields: 

y( X) = 2n(k) - n<**R t 1]-RC 1*1] ) - »««| l\\ 



+ 2n(Jc) - n(Jt+R[ zj-Rl " " »™J liil-Rf 3 ) 

♦ 2n<*> - n(k+R[ 3]-R[ 3+i - iM 3+1 -R J 

♦ 2n(k) - n(k+R( 41-R[ 4+H) - Mk+Rt **V R < 4 JJ 



(continued) 
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- 


n(k+R 


6+1 


-R[ 6]) 


+ 


2n(X) 


- 


a(k+R( 7 


l-Rf 7+1]) 


- 


n (k+R 


* 7+1 


-R[ 7]) 


+ 


2n(k) 


- 


a(k+R[ 8 


|-R[ a+i]) 


- 


n(k+R 


8+1 


-R[ 8]) 


+ 


2a(k> 


- 


a(k+R[ 9 


|-R[ 9+1]) 


- 


n(k+R 


; 9+1 


-R[ 9]) 


+ 


2a<k) 


- 


a<k+R{l(r 


|-R[10+1]) 


- 


a (k+R 


10+1 : 


-R[10]) 


+ 


2a(k> 


- 


a(k+R[ii 


|-R[11+1]) 


- 


a (k+R 


11+1 


-R[ll]) 


+ 


2n<k) 


- 


a(k+R[i2 


|-R[12+1]) 


• 


n(k+R 


12+1 : 


-R[12]) 


+ 


2a(k) 


- 


a(k+R(13 


|-R[13+1]) 


- 


a (k+R 


13+1 


-R[13]) 


+ 


2a<k) 


- 


a(k+R[l4 


|-R[14+1]) 


- 


a (k+R 


14+1 


-R[14]) 


+ 


2a(k) 


- 


a(k+R[15 


-R[15+l] ) 


- 


a (k+R 


15+1 ; 


-R[15]) 


+ 


2a(k) 


- 


a(k+R[16 


|-R[16+1]) 


- 


n (k+R i 


16+1 


-R[16]) 


+ 


2a(k) 




a(k+R[17 


|-R[17+1]) 




a (k+R 


17+1 


-R[17]) 


+ 


2a(k) 




a(k+R[18 


-R[18+l] ) 




a (k+R 


18+1 


-R[18]) 


+ 


2a(k) 




a(k+R[l9 


-R[19+l] ) 




a (k+R 


19+1; 


-R[19]) 


+ 


2a(k> 




a(k+R[20 


-R[20+l] ) 




a (k+R 


20+1 


-R[20]) 


+ 


2a(k) 




a(k+R[21 


-R[21+l]) 




n(k+Ri 


21+1 


•Rt2ll) 


+ 


2a(k> 




a(k+R[22 


-R[22+l] J 




a(k+Ri 


22+r 


-R[22]) 


+ 


2a(k) 




a(k+R(23 ; 


-R[23+l]) 




a (k+R 


'23+1 : 


-R[23]) 


+ 


2a(k) 




a(k+R[24 ; 


-R[24+l]) 




a(k+R| 


24+1 


-R[24]) 


+ 


2a<k) 




a(k+R[25 


-Rr25+1]) 




a(k+Ri 


25+1] 


-R[25]) 


+ 


2a(k) 




n(k+R(26; 


-R[26+l] ) 




a(k+R| 


26+1] 


-R[26]) 



(Eq A26) 

Evaluating R[j] in Equation (A26) using the definition of RJJ] found in Equation (A8) yields: 



y(k) = 


2D(k) 




a(k+l-2) 




a(k+2-l) 


+ 


2n(k) 




n(k+2-3) 




a(k+3-2) 


+ 


2a(k) 




n(k+3-0) 




a(k+0-3) 


+ 


2n(k) 




n(k+o-i) 




a(k+l-0) 


+ 


2&<k) 




a<k+l-2) 




a(k+2-i) 


+ 


2n<k) 




n(k+2-3) 




a(k+3-2) 


+ 


2n<k) 




n(k+3-0) 




a(k+o-3) 


+ 


2n(k) 




n(k+0-l) 




a(k+l-o) 


+ 


2a<k) 




a(k+i-2) 




a(k+2-l) 


+ 


2&<k) 




a{k+2-3) 




a(k+3-2) 


+ 


2a(k) 




a(k+3-0) 




a(k+0-3) 


+ 


2a(k) 




a(k+0-i) 




a(k+l-0) 


+ 


2n(k) 




a(k+l-2) 




a(k+2-l) 


+ 


2n(k) 




a<k+2-3) 




a(k+3-2) 


+ 


2a(k) 




afk+3-0) 




a(k+0-3) 


+ 


2n(k) 




a(k+o-i} 




a(k+l-0) 


+ 


2n(k) 




a<k+l-2) 




a(k+2-l) 


+ 


2n<k) 




a(k+2-3) 




a<k+3-2) 


+ 


2n(k) 




a(k+3-0) 




a(k+o-3) 


+ 


2n(k) 




a(k+o-i) 




a<k+l-0) 


+ 


2a<k) 




a<k+i-2) 




n(k+2-l) 


+ 


2n(k) 




a<k+2-3) 




a(k+3-2) 


+ 


2a<k> 




a(k+3-0) 




a(k+0-3) 


+ 


2B(k) 




a<k+o-i) 




a(k+l-0) 


+ 


2n<k) 




Q(k+l-2) 




a(k+2-l) 


+ 


2B(k) 




a(k+2-3) 




n(k+3-2) 



By combining similar terms, Equation (A27) can be reduced to: 

y(k) = 52n(k) - 20n(k-1) - 20n(k+1) - 6n(k + 3) - 6n(k-3) (Eq A28) 
Equation (A28) expresses the relationship between the measured data (as represented by the value of 
y(k) defined by Equation (A25)) and the estimated noise (as represented by the right side of Equation (A28)) 
as a function of the data frame parameter k. Since this relationship is true for ail the data collected, that is for 
all values of k, then the relationship of Equation (A28) can be rewritten in matrix form to cover all the data col- 
lected, as was done in Equation (C26): 

y = T • n , (Eq A29) 
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where £ and n are as expressed in Equations (C27) and (C28), respectively, and 



10 



15 



T = 



I 
i 

\ 



52 


-20 


0 


-6 


0 


0 


0 


20 


52 


-20 


0 


-6 


0 


0 


0 


-20 


52 


-20 


0 


-6 


0 


-6 


0 


-20 


52 


-20 


0 


-6 


0 


-6 


0 


-20 


52 


-20 


0 


0 


0 


-6 


0 


-20 


52 


-20 


0 


0 


0 


-6 


0 


-20 


52 



(Eq A31) 



Equation (A29) can be inverted to solve for the estimated noise matrix n; 
20 n = T- 1 ■ y , (Eq A33) 

where 

T _1 = the inverse of matrix T of Equation A31 . 

For example, if the selected window size w is 4, then the Toeplitz matrix T of Equation (A31) reduces to 
the (4 x 4) matrix: 

25 



30 



/ 
i 



T = 



52 


-20 


0 


-6 


20 


52 


-20 


0 


0 


-20 


52 


-20 


-6 


0 


-20 


52 



(Bq A3 6) 



35 The inverse T- 1 of this matrix T is; 

i 
I 



40 



1/252992 



6188 


3030 


1690 


1364 


3030 


7371 


3485 


1690 


1690 


3485 


7371 


3030 


1364 


1690 


3030 


6188 



(Eq A371 



45 Thus, for example, the (k^k^ = (2,3) term of the matrix T- 1 : 

T-KMz) = T-1(2,3) = 3485 / 252992. (Eq A38) 
For the filter window size w = 4, Equation (C31) can be expressed as: 



50 * < 

n<X> = 1/4 Z Z (T 'Oc,,*,) . y(k+K T -)c 2 J ), 



(Eq A38a) 



where the matrix elements T- 1 (k 1t k2) are given by Equation (A37) and the terms y(k+k r k 2 ) are given by Equa- 
55 tion (A25). 
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Minimal Discontinuity Filter for Alternative Sampling Rate 

As mentioned earlier, the minimal discontinuity filter derived in the previous section was based on a sam- 
pling interval d of 0.4 iches. If, for example, the sampling interval d is 0.1 inches instead, so as to give better 
vertical image definition, then the general function g(j), which expresses the relationship between the data 
frame parameter k and the hole depth parameter i and the electrode parameter j, is: 

k = i - g(j) = i - 4R[j], (Eq A39) 
where Rfl is defined by Equation (A8). 

The filter derivation then follows a path similar to that for the 0.4-inch sampling depth interval, except that 
wherever the functions R[j] and R[j+1] appear in Equations (A7) ( (A24), (A25), and (A26), the functions 4R[j] 
and 4R[j+1] should be substituted, respectively. The result is the following equation relating the measured data 
in the form of y(k) to the estimated noise n(k), which is analogous to Equation (A28): 

y(k) = 52n(k) - 20n(k-4) - 20n(k + 4) - 6n(k+12) - 6n(k-12) (Eq A40) 

where 



15 



20 



26 



M<k+4R[j] , j+l) 



y(k)= t <K<*MR[j],j> 

j=l 

- K<K+4R[j+l] , j) + H<*+4R[j+l] , j+l> > 



(Zq A41) 



25 



This relationship can still be presented in the matrix form of Equation (A29): 

y = Tn, (Eq A42) 
where the Toeplitz matrix T still has the form expressed in Equation (A31), but now: 



30 



35 



40 



1 = 



/ : \ 

YU-12) 

Y<*-8) 
Y<*-4) 

Y<*) 
Y<*+4) 
Y<*+8) 

Y<*+12) j 

i 

\ : / 



(Eq A43) 



45 



50 



and 
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10 



15 



20 



25 



{Eq A44) 



/ : 
i 

n<k-i2) 
n(X-8) 
n(X-4) 
n = : n<X) 
n(X+4) 
nU+8) 
n<k+l2) 



The minimal discontinuity filter for this example has a form analogous to that expressed in Equation (C31): 



w v 

n<)c> = 1/w z r [T' l (X x ,k 7 ) 



(Eq A45) 



Thus, for a given value of k (that is, for a given data frame), only measured data from every fourth data 
frame need be included in the computations for the estimated noise n(k) of that data frame k. This makes sense 
when it is remembered that for an array with a = 0.4 inches between electrode rows and a sampling interval 
d of 0.1 inches (as in this example), the electrodes of Row 2, for instance, will reach the hole depth of the elec- 
30 trodes of Row 1 only after 4 data frames have been taken. In essence, the current example is one in which 4 
linearly independent filters offset from each other by 0.1 inches are essentially operating in parallel. 

Application of the Minimal Discontinuity Filter 

35 The general equation for the minimal discontinuity filter of this invention used to approximate the one-di- 
mensional noise n(k) is given by Equation (C31) as: 



40 



n<k> = 1/w Z z (T* l (k lf k t ) 



<Eq Bl) 



45 



50 



where the function y(k) has the general form as given in Equation (C25): 



E-l 

Y<X) = £ <K(X+g<j),j) - M(k+q(j) , - M(X+q(j+l) , j) + 
j+1 

M(k+g<j+l) , } (Eq B2) 



where M(k+g(j),j) = M(i,j) = the data measured by electrode j at depth i during data frame k. 

One particular array that can be used in the logging of bore holes is the one described in the previous sec- 
tions and shown in Figure 6. The implementation of the minimal discontinuity filter invention can be demon- 
55 strated using the set of hypothetical data shown in Figure 7 and assumed to be collected during a bore-hole 
pull using the array of Figure 6, but wherein data is collected every 0.4 inches. As can be seen from this hy- 
pothetical data, when Row 1 of the array was at a depth of 8.4 inches, the one-dimensional noise occurred 
such that an additional 20 units were added to each of the 27 electrode signals. For example, the (j=18) elec- 
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trade in Row 3 should have measured 6 units to correspond to the true characteristic of the hole at that depth, 
but because of the linear noise of 20 units, the recorded signal was 26 units. 

This example will demonstrate the implementation of the minimal discontinuity filter of Equation (B1) by 
computing the noise contained in the signal measured by electrode j=18 at a depth of 9.2 inches. The parameter 
5 k in Equation(BI) is used to identify the frame of 27 electrode signals measured at the same time; the para- 
meter i is used to identify the set of 27 electrode signals corresponding to the same depth in the hole. In the 
present example, the depth of 9.2 inches is identified by the parameter i=23. 

For the array of Figure 6 and a sampling interval of 0.4 inches, the relationship between k and i and j as 
expressed by the function g(j) is given in Equation (A7) by: 
10 k = i - g(j) = i - R[j], (Eq B3) 

where R[j] is as defined in Equation (A8). Thus, for electrode j=18 at depth i=23, the data frame parameter k 
is found using Equation (B3): 

k = i - Rfj] = 23 - R[18] = 23 - 2 = 21, 
the depth of the top row of the sawtooth artifact depicted in Figure 7. 
15 For the array depicted in Figure 6 and a sampling depth interval of 0.4 inches, the Toeplitz matrix T is given 
by Equation (A31). For this particular calculation, a filter window size w of 4 is used. When the filter size is 4, 
the Toeplitz matrix of Equation (A31) reduces to the (4 x 4) matrix of Equation (A36). The inverse T- 1 of this 
matrix T is given by Equation (A37). 

The one-dimensional noise in electrode j=18 at depth i=23 is estimated by computing n(k) = n(21) using 
20 Equation(BI) (where the filter window size w for this example is set to 4): 



4 4 

n(k) =1/4 1 E (T'M**,**) yOc+Xj-XiU (Eq B8) 

25 *,=1 )Cj=l 

Equation (B8) can be expanded to yield: 

30 

n(k) = 1/4 

(T' l (l,i) y(k) + T" x (l,2) Y()C-1) + T l (l,3> y(X-2) 

+ T' l C2,l) y<X+l) + T* 1 (2 / 2> y()c) + T'^Z, 3) y(*-l) 

+ T -1 {3,1) y<lt+2) + T -i (3,2) y<X+l> + r l <3,3> y(K) 

+ T" l (4,l) y<*+3> + T- x (4,2) Y<*+2) + T" l <4,3> y(*+l) 

35 



Substituting 21 for k yields: 

40 n(X) = n(2l) = 1/4 

<T' l (l,l> y(2l) + T^d^) y(20) + T' l (l,3) y(19) 
+ T* l (2,l) y<22) + T* i (2 # 2) y(21) + T* A <2,3) y(20) 
+ T* l (3,l) y(23) + T' l (Z,Z) y<22) + T _l <3,3) y(21) 
+ T' l (4,l) y(24) + T' A (4 # 2) y(23J + T" l (4,3) y<22) 

45 

For the array of Figure 6, the function y(k) of the measured data in Equation (B1) is found by substituting 
R[j] for the function g(j), as expressed in Equation (B3), into Equation (B2), and by noting that E-1=26 for an 
array with 27 electrodes, to yield: 

50 

26 

y(k) = I {K(k*R[j],j) - H<*+R[j], j+1) + M(Jc+R[j+l] ,j+l) - 
j=l 

M(k+R[j + 1], j) > (Eq BID 

55 

Equation (B11) can be expanded to yield: 



+ T _1 (l,4) y(Jc-3) 
+ T"M2 f 4) y(lt-2) 
+ T" l (3,4) y(k-l) 
+ T' l (4,4) Y<*)> 

(Eq B9) 



+ T' l (l,4) y(18) 
+ T' l (2,4) Y(19) 
+ T* l (3,4) y(20) 
+ T* l (4,4) y(2D) 

(Eq BIO) 
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= 


H(k+1, 


I) 


• 


M(k+1, 


2) 




M(k+2, 2) 


— 


MCfc+2, 


1) 


+ 


M(X+2 , 


2) 




M(k+2, 


3) 


+ 


M(X+3, 3) 




M(X43, 


2) 


+ 


K(X+3, 


3) 


• 


M(k+3, 


4) 




M(X , 4) 




M(X , 


3) 


+ 


M(X , 


4) 


— 


M<k , 


5) 


+ 


M(X+1, 5) 




M<X+1, 


4) 


+ 


M(X+1, 


5) 


— 


M<k+1, 


6) 




M<X+2, 6) 


** 


M(X+2, 


5) 


+ 


H(X+2, 


6) 




K(k«, 


7) 


+ 


M(X+3, 7) 




M(X+3, 


6) 




H(X+3, 


7) 




M(k+3, 


8) 




M(X , 8) 




M(X , 


7) 


+ 


M(X , 


8) 




K(k , 


9) 


+ 


K(X+1, 9) 




M(X+1, 


3) 


+ 


M(X+1, 


9) 




K(k+1, 


10) 


+ 


K(X+2,10) 




M<X+2, 


9) 


+ 


M0C+2, 


10) 




H(k+2, 


11) 




K(X+3,11) 




M(X+3, 


10) 


+ 


M(X+3, 


11) 




M(k+3, 


12) 


+ 


M(X ,12) 




M(X , 


11) 


+ 


M(X ,12) 




M(k , 


13) 




KOt+1,13) 




M(X+1,12) 


+ 


M(k+1, 


13) 




H(k+l,14) 




H(X+2,14) 




M(X+2,13) 


+ 


M(k+2, 


14) 




M(k+2, 


15) 




M(X+3,15) 




MOC+3, 


14) 


+ 


M(k+3, 


15) 




M(k+3, 


16) 




H(X ,16) 




M(X ,15) 



(continued) 



+ M(k ,16) 
♦ H(X+i,i7) 
+ M(X+2,18) 
+ M(k+3,19) 
+ M(X ,20) 
+ M{X+1,21) 
+ M{X+2,22) 
t M(X+3,23) 
+ M(X ,24) 
+ M(X+1,25) 
+ M(X+2,26) 



- M<X ,17) 

- K(k+l,18) 

- K(X+2,19) 

- K(X+3,20) 

- M(k ,21) 

- K(X+1,22) 

- M(X+2,23) 

- M(X+3,24) 

- M(X ,25) 

- M(X+1,26) 

- M(X+2,27) 



+ M(X+1,17) 
+ M(X+2,18) 
+ M(X+3,19) 
+ H<X ,20) 
+ M(X+1,21) 
+ M(X+2,22) 
+ M(X+3,23) 
+ M(X ,24) 
+ K(X+1,25) 
+ H<X+2,26) 
+ M(X+3,27) 



- MCX+1,16) 

- K(X+2,17) 

- K(X+3,18) 

- M(k ,19) 

- K(X+1,20) 

- K(X+2,21) 

- M(X+3,22) 

- M(X ,23) 

- H(X+1,24) 

- M(k+2,25) 

- H(X+3,26) 



The expansion of Equation (B1 0) requires the computation of y(k) using the expansion of Equation (B12) 
seven times, for values ofk= 18, 19, 20,21,22, 23, and 24. Fork= 18, Equation (B12) becomes: 



y(i8) 



= M(19, 
+ M(20, 
+ M(21, 

+ M(ia, 

+ KU9, 
+ M(20, 
+ M<21, 
+ M(18, 
+ M(19, 
+ M(20. 



1) 
2) 
3) 
4) 
5) 
S) 
7) 
8) 
9) 
10) 



+ M<21,11] 

+ M(ie r i2) 

13) 
14) 
15) 
16) 
17) 

18) 

+ M<21,19) 
+ M<18,20) 
+ M(19,21) 
+ M(20,22) 
+ 8(21,23) 
+ M(18,24) 
+ M(19,25) 
+ M(20,26) 



+ M(19 
+ M<20 
+ M(21 
+ M(18 
+ M(19 
+ M420 



4) 

5) 
6) 
7) 



H(19, 2) 
M(20, 3) 
M(21, 
M(18, 
H(19, 
K<20, 
M{21, 8) 
H(18, 9) 
M(19,10) 
K(20,ll) 
M(21,12) 
K<18,13) 
M(19,14) 
M<20,15) 
K<21,16) 
M(18,17) 
M<19,18) 
M(20,19) 
M(21,20) 
M(18,21) 
M<19,22) 
M(20,23) 
M<21,24) 
M(16,25) 
M(19,26) 
M(20,27) 



2) 
3) 
4) 
5) 
6) 
7) 
8) 
9) 



K(20, 
M(21, 
M<18, 
M(19, 
H<20, 
M(21, 
K(18, 
M(19, 
M(20,10) 
M(21,ll) 
M(18,12) 
M(19,13) 
M(20,14) 
M(21,15) 
M(18,16) 
M(19,17) 
K<20,18) 
M<21,19) 
H<18,20) 
M(19,21) 
M(20,22) 
M(21,23) 
M(18,24) 
M(19,25) 
M|20,26) 
H(21,27) 



1) 
2) 
3) 
4) 
5) 
$) 
7) 
8) 
9) 



H<20, 
M(21, 
M(18, 
M(19, 
M(20, 
M<21, 
M(18, 
M(19, 
M(20, 
M(21,10) 
M(18,ll) 
M(19,12) 
M(20,13) 
M(21,14) 
M(18,15) 
H(19,16) 
M(20,17) 
M(21,18) 
M(1B,19) 
M(19,20) 
M(20,21) 
M(21,22) 
M(18,23) 
M(19,24) 
M(20,25) 
M(21,26) 



(Eq B13) 
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The values of this expansion are found in Figure 7. For example, the measured signal M(19,1) is the signal 
measured by electrode j=1 at depth i=19, or 2 units. The signal M(19,2) measured by electrode j=2 at depth 
i=19 is also 2 units. Note, for example, that the first two terms of the third row of the expansion have values 
M(21,3) = 4 units and M(21,4) = 24 units. Substituting the appropriate values from Figure 7 into Equation (B13) 
5 yields: 



y(i8) 



10 



15 



20 



25 



30 



= (2 

+ 3 

+ 4 

+ 1 

<► 2 

+ 3 

+ 4 

+ 1 

+ 2 

+ 3 

+ 4 

+ 1 

+ 2 

+ 3 

+ 4 

+ 1 

+ 2 

+ 3 

+ 4 

+ 1 

+ 2 

+ 3 

+ 4 

+ 1 

+ 2 

+ 3 



2 
3 
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Similar evaluations of the expansion of Equation (B12) can be made to solve for the other 6 y(k) values in Equa- 
tion (B10) to yield: 



y(l9) s o units 

y(20) = -400 units 

y<21) = 1040 units 

y<22) = -400 units (Bq B15) 

y(23) = o units 

y{24) = -120 units 



Substituting these values fory(k) together with the values for the sensors of the inverse matrix T- 1 of Equation 
45 (A37) into the expansion of Equation (B10) yields: 



n(fc) = n(21) = (1/4) (1/252992) 

{ («188) (1040) + (3030) (-400) 
+ (3030) (-400) + (7371) (1040) 
+ (1690) (0) + (3485) (-400) 
+ (1364) (-120) + (1690) (0) 

=20 units. 



* (1690) (0) + (1364) (-120) 

+ (34B5)(-400) + (1690) (0) 

+ (7371) (1040) + (3030) (-400) } 

+ (3030) (-400) + (6188X1040)} 



This computation shows that the noise n(k) when k is 21 is estimated to be 20 units. The measured signal M(i,j) 
55 is the sum of the ideal, true signal S(i,j) (corresponding to the true value of the measured characteristic of the 
hole at depth i and electrode j) and the noise N(ij) (corresponding to the noise at depth i and electrode j), or 

M(ij) = S(i t j) + N(i,j) (Eq B16) 
The true signal S(iJ) is therefore expressed by: 
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s(ij) = M(i,j) - N(ij) = M(ij) - n[i-g(j)] = M(ij) - n(k) (Eq B17) 
The one-dimensional noise which this minimal discontinuity filter invention addresses is such that all 27 
electrodes are affected at the same time by the same amount. An estimated noise value n(k) of 20 units for 
k=21 means that all electrodes with k value 21 should have 20 units substracted from their measured signals 
5 M(ij) in order to yield the true data signal values S(i,j). Equation (B3) can be used to identify the electrode 
signals which correspond to k=21. They are: 

M(22,1), M(23,2),M(24,3), M(21 ( 4) M(18 f 23)„.., M(23,26), M(24,27) 

The estimation of the true signal S(ij) for electrode j=18 at depth i=23 using Equation (B17) is therefore: 
S(23,18) = M(23 ( 18) - N(23,18) = M(23,18) - n(21) = 26 - 20 = 6 units. 
10 In this example, using ideal hypothetical values, the minimal discontinuity filter of Equation (B1) produced a 
value exactly equal to the amount of assumed noise which in this case is 20 units, thereby yielding a value for 
the estimated true signal exactly equal to the true value of the measured characteristic of the hole. 

The computations in this example were geared towards solving for the noise level in one particular elec- 
trode sampling. Under normal circumstances, the flow of mathematical processing is preferably slightly dif- 
15 ferent A more realistic example of the use of the array of Figure 6 for bore-hole logging may involve as many 
as 10,000 frames of data and a filter window size of 11. The sequence of computations on the measured data 
would then be: 

(1) Compute the inverse T- 1 of the (11 x 11) Toeplitz matrix T having the form of the matrix in Equation 
(B2), 

20 (2) Compute the function y(k) of the measured data using Equation (B11) for data frames k = 1 to 10,000, 

(3) Compute the estimated noise level n(k) using Equation (B1) for data frames k = 1 to 1 0,000 and filter 
window size w = 11 , and 

(4) Calculate the estimated true electrode signals S(ij) for electrode j and depth i using Equation(B17) for 
j = 1 to27andi = 1 to 10,000. 

25 Figures 8, 9, and 1 0 show, respectively, a flow diagram for the broad steps of the method of the invention, 
a flow diagram for the steps involved in one preferred form of the method, and a more detailed flow diagram 
for a preferred method according to the invention, the implementation of which may be accomplished by ap- 
propriate corresponding programming of the minicomputer 54. 

Figure 8 shows at 100 that one first defines the g(j) of the selected array, as described above. At 102 the 

30 depth-adjusted sample values M(i,j) are accessed, together with their identifiers i and j. At 104 there is gen- 
erated the correcting signal value of n[i-g(j)] which minimizes a measure I of the image discontinuity which is 
a function of differences transverse to the logging direction: S(ij) - S(i,j+1), or M(j) - n[i-g(j)] - M(i,j+1) + n[i- 
gfl+1)]. At 106 each value of the correcting signal n[i-g(j)] is subtracted from M(ij) to produce the corrected 
values S(i,j). At 108, the corrected signals S(iJ) are applied to the display device to form the two-dimensional 

35 image. 

In the somewhat more detailed flow diagram of Figure 9, at 11 0 and 112 the first two steps of Figure 8 are 
repeated, and at 116 is shown the use of M(i,j) and g(j) to form the set of simultaneous equations resulting 
from taking the partial derivative of I with respect to n. At 118 is shown the formation of the matrix equation T 
• n = y, corresponding to d\l 5n[i-g(j)]. At 120 the inverse T _1 of matrix T is formed, for a window w, selected at 
40 122. At 124 the correcting signal values n(k) are formed from J- 1 and y. The data is corrected and then used 
to produce the two-dimensional image as already shown at 106 and 108 of Figure 8. 

Figure 10 shows the preferred process in detail. In this case, since the expression for n(k) shown partly 
at 130 and partly at 134, is known from the teachings of this invention, it need not be derived again in practical 
applications of the method; all that is necessary is to substitute into this expression the relevant values of g(j) 
45 and E to evaluate it for a particular case. This substitution is shown at 130 and 132, respectively. One first 
defines the form of g(j) at 134 for the particular array and sampling interval; for example, for the array of Figure 
6, g(j) this may be 0.4 R[j], using a sampling depth interval d equal to 0.1 inch and an inter-row distance of a 
= 0.4 inch. The total number of electrodes E is defined at 136, which number in the foregoing example is 27. 

These values of g(j) and E are substituted into the portion of the expression for n(k) which contain no meas- 
50 ured signal values M(i,j), as shown at 130. They are also substituted into the expression y for the terms con- 
taining M(i,j), as shown at 132. 

Considering first the steps flowing from 130, the expression in 130 is expanded for j = 1 to j = (E-1) at 138. 
The result is evaluated for g(j) and g(j+1), at 140. At 142, terms are combined and simplified, and the results 
used to generate the matrix T, at 144. Next the window value w is selected at 146 and used to define the (w 
55 x w) curtailed matrix T. The latter matrix T is used to generate its inverse T~ 1 at 150, and T~ 1 is then used in 
the generating of n(k) at 152. Once T- 1 has been computed for the array and the system to be used, it need 
not be computed again for that array and system. 

Steps flowing from 132 involve substitution, into the formula shown, of the defined g(j) and the defined 
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E; the y(k) so generated are then expanded for j = 1 to (E-1), at 160. At 162, D is defined equal to the number 
of data frames to be examined (e.g. 10,000). At 164, M(i,j) are accessed, and at 166 the expanded y(k) use D 
and M(ij) from 162 and 164, respectively, to generate y(k) for j = 1 to D. 

At 1 70, the selected value of w is substituted into the expression for n(k) and at 1 72 the result is expanded 
5 for ki = (1 to w) and for k 2 = (1 to w). At 176 the expanded value of n(k) is used to compute n(k) for k = (1 to 
D), using T- 1 from 150 and y(k) from 166. 

The resultant n(k) from 176 and the M(i,j) signals from 164 are used at 178 to correct the M(ij) values by 
subtracting n(k) from M(i,j) for i = 1 to D and j = 1 to E. As indicated at 180, the resultant corrected signals S(i,j) 
are used in forming the final two-dimensional image. 
10 Figure 11 shows the minimal discontinuity filter in apparatus form. The information g(j) and E is supplied 
to a matrix generator 200 which generates the matrix T described above. From T, the matrix inverter 202 pro- 
duces T" 1 as defined previously, and supplies it to the noise estimator 204. The latter estimator is also supplied 
with the g(j) and E information, as well as the measured data M(i,j) and D, and from these produces the esti- 
mated noise signals n(k) as also previously described. The signal corrector 206 subtracts n(k) appropriately 
15 from M(i,j) to produce the corrected S(i,j), which are used to generate, in the image generator 208, the signals 
which produce the actual two-dimensional bore-hole image on the image-display device. 

Fig. 12 shows at A and B two logs taken on the same logging run by two different pads of the formation 
microscanner referred to above, without filtering according to the invention. Fig. 12 shows at C and D respec- 
tively the logs obtained by applying the filtering system of the invention to the same data which produced logs 
20 A and B, and in which the sawtooth component is greatly reduced or eliminated. 

In the idealized example discussed above with respect to Figure 7, the true signal values in each horizontal 
row are assumed all to be equal, and the noise exists for exactly one frame time and is of exactly constant 
value. Under these conditions the correcting signal produced exact and complete cancellation of the assumed 
noise. To the extent that these ideal conditions are not present in a particular situation, noise cancellation will 
25 not necessarily be exact, and the measured values after noise cancellation of samples obtained when no noise 
is present may differ somewhat from their true values; nevertheless, for a substantial range of departures from 
ideal conditions, the correcting signal will provide useful reduction of the artifacts described above. 

There are a variety of other respects in which the method and apparatus of the invention may differ sub- 
stantially from that illustrated. For example, the noise may start and/or end in mid-frame; or there may be a 
30 plurality, or a large number, of immediately successive frames during which the noise occurs; or there may be 
a large number of both closely and widely time-spaced periods of noise, and in all such cases the noise may 
well differ in intensity during different frames. It is also assumed in the foregoing examples that the system is 
linear in the sense mentioned above, so that M = N + S, e.g., it is not so heavily overloaded by N and/or S that 
substantial non-linearity occurs. 
35 With regard to the sensors, they need not be electrodes for an FMS™ or SP system, but can constitute 
any means for deriving values indicative of a characteristic of the adjacent portions of the earth formation, using 
for example acoustics, magnetics, particle emission, etc. 

The pattern of the sensors in the array can also vary greatly in different embodiments of the invention. 
Any number of sensors greater than one may be used, and they may be in any pattern so long as at least some 
40 of them are spaced from each other along the direction of array motion and not aligned with each other along 
that direction. The array can be described by sets of sensors, where each set is identified by a position along 
the direction of array motion. 

For example, one possible array uses a plurality of sensors arranged along a straight line not parallel to 
the direction of array motion, e.g., one ramp of the previously-described sawtooth. In such an array, the number 
45 of sets of sensors equals the total number of sensors in the array, each set containing just one sensor. 

In another situation, the logging probe may contain sensors which are not all in unique locations along a 
second direction transverse to the direction of array motion; that is, there may be at least two sensors aligned 
along the direction of array motion. Under such conditions, the sets of sensors which comprise the array for 
purposes of this invention are those which are not aligned along the logging direction. Thus, for example, a 
50 probe which contains four sensors arranged in a square pattern with two of the sides of the square parallel to 
the logging direction will have two pairs of sensors which are aligned along that logging direction. In such a 
case, the array of sensors for the purpose of this invention may be defined by two sets of sensors, the first 
set containing only one of the sensors at the top of the square pattern and the second set containing only the 
sensor at the bottom of the square which is not aligned with the sensor of the first set, that is, the sensor at 
55 the other end of the diagonal bisector of the square pattern. This definition of the array of sensors includes 
two of the sensors of the probe and excludes the other two, for the purposes of this invention. 

The parameters i, j, and k need not be used to identify the array position and the positions of the sensors 
according to the system described above, so long as k in some way identifies the array position along the log- 
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ging direction and i and j identify the electrode positions according to a g(j) which permits expression of i in 
terms of j. Nor is it necessary to measure i and k in terms of increasing depths, as opposed to decreasing 
depths. 

With regard to the pulling rates, sampling intervals, frame durations, and spacings of sensors along the 
5 pulling direction, the system is greatly simplified by performing the sampling so that the array does not move 
appreciably during the time it takes to complete one sampling, and by using the "frame grabber" illustrated, 
all sensors can in fact be sampled at exactly the same time; if instead the sensors themselves are sampled 
sequentially, each sampling frame is preferably of very short duration, i.e., the sensors are scanned at a very 
high rate. Further, while the pulling rate, the sampling time interval, and the distances between sensors along 
10 the pulling direction may be coordinated so that all sensors tend to be sampled for each depth at which any 
samples are taken, as pointed out above this is not essential, nor is it absolutely essential that the spacings 
between sensors along the pulling direction all be the same, since the signal processor can accommodate a 
wide range of variations in these respects. 

It is also noted that the filter is not limited to use with data from the logging of a bore-hole wall in the pro- 
fs specting for oil; the wall can be of any shape, size, and origin, and explored in searching for different minerals 
or for scientific geologic purposes, as examples. Also, the logging motion need not be exactly vertical, although 
the rows of sensors are preferably at right angles to whatever the direction of motion is. 

With regard to the filter of the invention, it will be understood that it will usually not be used in real time in 
the logging process; the logging data are typically obtained at one time, and may be examined and processed 
20 at leisure, at a later time. At such time, the raw data, the depth-adjusted data, and the data passed through 
the filter of the invention may ail be separately stored, and retrieved when desired for analysis and such proc- 
essing as is desired. In general, one skilled in the art will be able to recognize the presence of one-dimensional 
noise in the uncorrected image which will benefit from use of the minimal discontinuity filter, and wiit use the 
filter to reduce the noise- induced artifacts. In some cases, as where the noise persists for at least several 
25 frames, repeated application of the filter of the invention may be made, by recirculating the corrected samples 
so they are filtered one or more additional times. The window size, and the number of frames filtered, can 
also be varied by the user to obtain best and/or simplest filtering. 

While preferable in most cases, a window need not be used, and instead direct evaluation of a complete 
selected set of data may be used, especially where a high-powered computer is available. Where a window is 
30 used, it is preferably used only over a range of samples not so near the top or bottom of the surveyed area 
that a complete sliding of the window cannot be performed. 

Also, while the above-described specific mathematical expressions for the filter are preferably employed 
in producing the correcting signal n, variations in this procedure are possible so long as the filter minimizes 
discontinuities between samples corresponding to the same depth, and, in effect, senses the presence of one- 
35 dimensional noise in the depth-adjusted samples and reduces that noise by subtracting estimated values of 
noise from all samples, the subtracted value being the same for all samples produced in the same original frame 
of samples. 

For example, in the estimate of n(k) set forth above, the sum of the squares of differences S(i,j) - S(ij+1) 
was used as a measure I of the absolute value of differences in true signal values along the transverse direc- 
40 tion. However, other measures I may be used including, without limitation, the sum of higher even powers of 
S(i,j)-S(ij + 1). 

Further, while the filter of the invention finds special utility in use with signal samples from a staggered 
array, it is also useful in reducing certain types of noise which produce spurious differences between samples 
corresponding to positions on the bore-hole wall which are aligned with each other along the pulling direction, 
45 rather than staggered. 



Claims 

so 1 . A method for producing a two-dimensional image of a characteristic of an area of a wall (1 4) of an earth 
formation (16) using an array (12) of sensors (30) comprising a plurality of sets of one or more sensors 
(30), some of said sets being at different positions in said array (12) as measured along said first direction, 
and at least one of the sensors (30) of at least one of said sets being differently positioned from all of the 
sensors in at least one of the other of said sets, as measured along a second direction transverse to said 

55 first direction, the method comprising moving the array (12) along said wall (14) in a first direction to pro- 

duce a signal level at each sensor (30) representative of the value of said characteristic in the portion of 
said wall (14) adjacent thereto, sampling (36) said signal levels at all of said sensors (30) of said sets at 
each of a succession of sampling times, to produce successive original frames of sensor signal samples, 



25 



EP0 426 564 B1 



processing (60) said samples and applying said samples to a two-dimensional display device (68) in sets 
such that said processed samples of each of said original frames control points in said display which are 
in a geometric pattern corresponding to the geometric pattern of said sets of sensors (30) in said array 
(12), characterized by prior to applying said processed samples to said image-display (68), generating 
(69) from said samples a correcting signal n, and subtracting it from each of said samples, to reduce any 
artifact in said image produced by noise affecting said sensors, said noise affecting each of said sensors 
(30) substantially equally during any of said sampling times and producing an artifact in said image having 
a form corresponding to said pattern of said array (12) of sensors (30), wherein the values of n are gen- 
erated such as to minimize abrupt discontinuities between values of those samples which are produced 
in different original frames by sensors (30) at the same position along said first direction, and wherein 
the values of n so subtracted are equal for all samples produced at the same sampling time. 

A method as claimed in claim 1, wherein said correcting signal n has values which minimize an image 
discontinuity measure I, where I is a function of differences in signal samples evaluated over a plurality 
of said frames of samples. 

A method as claimed in claim 1 or 2, wherein said sensors (J1-J27) of said sets which are located at dif- 
ferent positions along said first direction are spaced apart along said first direction by distances (a) which 
are integral multiples, including one, of the distance (d) which said array (12) moves between successive 
sampling intervals. 

A method as claimed in claim 2, wherein I is a function of S(i, j) - S(i, j+1), where S(i, j) = M(i, j) - N(i, j), 
where: 

i is an integer index representing the position (a) of each of said sensors (30), with respect to a 
reference position, when its signal level is sampled, measured along said first direction and in terms of 
units equal to the distance (d) which said array moves between successive samplings; 

j is an integer index representing the position of each of said sensors (30) along a second direction 
transverse to said first direction, and measured in terms of the ordinal number of each sensor position, 
numbered along said second direction; 

M(i, j) is the measured value of each of said sensor signal samples taken at sensor positions (i, j); 

N(i, j) is the value of said noise in the signal sample at each corresponding sensor position (i, j) 
and is substantially the same for all of said sensors (30) at any given sampling time; and 

S(i, j) is the true value of the sensor signal sample produced by each of said sensors at position 
(i, j), in the absence of said noise N(i,j). 

A method as claimed in claim 4, wherein said image discontinuity function I is a function of the absolute 
values of said differences S(i, j) - S(i, j+1), evaluated over a plurality of said frames of samples. 

A method as claimed in claim 4 or 5, wherein said image discontinuity function I is a function of the sum 
of the squares of said differences S(i, j) - S(i, j+1), evaluated over a plurality of said frames of samples. 

A method as claimed in claims 4-6, wherein said correcting signal n has values satisfying the set of si- 
multaneous linear equations derivable from d,/d n = 0. 

A method as claimed in claim 7, wherein said correcting signal has values equal to those resulting from 
an evaluation of. 



D 

n(k)= I T" 1 (k.x).y(x), 
x=l 

wherein: 

n(k) is the estimated value of said noise in each of said sensor signal samples measured during 
the same said frame k of sampling, 

x is an integer summation parameter which takes on values from 1 to D, 
D is the number of frames comprising substantially the entire image, 
"T 1 (k. x) represents the (k, x) element of the matrix T~ 1 , 
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T _1 represents the inverse of a banded symmetric Toeplitz matrix T, 

T represents in matrix form the coefficients of the noise-containing terms in the set of linear si- 
multaneous equations obtainable by setting d\Id n = 0, 



y(x) = 1 (M(x+g(j), j) - M(x+g(j), j+1) - M(x+g(j+l). j) + M(x+gG+l), j+DJ 
J =1 

10 E represents the number of said sensors (30) in said sensor array (12), and 

g(j) represents in units of i the position of each of said sensors (30) along said first direction with 
respect to a reference position fixed with respect to said array (30), expressed as a function of the position 
of each sensor (30) along said second direction (j). 

15 9. A method as claimed in claim 7, wherein said correcting signal has values equal to those resulting from 
an evaluation of: 



w w 

n(k) = 1/w I I T'Hkhkl) • y (k+ki - k2) 
20 ki=l k2=l 

wherein: 

n(k) is the estimated value of said noise in each of said sensor signal samples measured during 
25 the same said frame k of sampling, 

w represents the size of the filter window which is equivalent to the dimension along said first di- 
rection of the size of said portion of said image, where w is some number less than the number of frames 
which comprise the entirety of the image, 

kj and k 2 represent the integer summation parameters which take on values from 1 to w, 
3 o T-^k^ kj represents the (k 1p kj element of the matrix T~ 1 , 

T-1 represents the inverse of a banded symmetric Toeplitz matrix T, 

T represents in matrix form the coefficients of the noise-containing terms in the set of linear si- 
multaneous equations obtainable by setting djd n = 0, 

y(k) represents the function of said measured data, whereby: 

35 

E-l 

y(x) = I {M(x+g(j), j) - M(x+g(j), j+l) - M(x+g(j+l), j) + M(x+g(j+l) f j+l)J 

40 

E represents the number of said sensors in said sensor array, and 

g(j) represents in units of i the position of each of said sensors (30) along said first direction with 
respect to a reference position fixed with respect to said array (30), expressed as a function of the position 
of each sensor (30) along said second direction (j). 

45 

10. A method as claimed in claim 9, wherein said integral number is 4 and said function g(j) is MRQ], where 
RQ] is the remainder resulting from division of j by 4. 

11. A method as claimed in any preceding claim, wherein said signal samples are depth-adjusted before sub- 
traction of said correcting signal n by selecting, from the frames of samples produced by sensors at dif- 
ferent depth positions, corresponding sets of samples produced by said sensors when they are at the 
same depths. 

12. A system for producing a two-dimensional image of at least a segment of a surface of an earth-formation 
wall (14), comprising: 

55 (a) a two-dimensional array (12) of sensor electrodes (30) and means (27,26) for moving said array 

(12) along said wall segment (14) in a first direction, the center of said electrodes (30) being spaced 
from each other both along said first direction and along a second direction at right angles to said first 
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direction, at least some of said electrodes (30) at different positions along said first direction being non- 
aligned with each other along said first direction; 

(b) sampling means (36, 56) for repetitively sampling the signal levels at ail of said electrodes (30) to 
form corresponding original frames of signal samples; 
5 (c) depth-adjustment means (60) for processing said original frames of signal samples to produce 

depth-adjusted frames of said signal samples, all of said samples in each of said depth-adjusted frames 
being produced at the same position along said wall (14); 

d) a two-dimensional display device (60) responsive to said depth-adjusted frames of signal samples 
to produce successive lines in said image in response to corresponding different ones of said depth- 
10 adjusted frames of signal samples, thereby to form a two-dimensional image of said wail segment (14), 

characterised by: 

(e) filter means (62) for filtering said depth-adjusted signal samples to reduce artifacts introduced 
into said image by noise which affects all of said electrode signal levels substantially equally, said artifacts 
corresponding to the patterns of electrodes (30) in the array (12), said filter means (62) comprising noise 

15 signal estimating means responsive to said depth-adjusted samples to generate correcting signals having 

discrete values n(k) identified with each of said depth-adjusted samples, and means for subtracting each 
value n(k) of said correcting signals from the corresponding sample to reduce said artifacts in said image, 
and wherein said correcting signal values are selected to minimize differences between the values of said 
depth-adjusted signal samples produced by electrodes adjacent to each other along said second direc- 

20 tion, and to distribute said correcting signals equally over sets of depth-adjusted samples produced by 

the same original frames of said signal samples. 

13. A system as claimed in claim 12, wherein said values N(k) are generated so as to minimize the differences 
S(i, j) - S(i, j+1), where S represents the true values of said depth adjusted signal samples in the absence 

25 of said noise, I is an index representing the depths at which said depth-adjusted signal samples were orig- 

inally produced, and j is an index representing the ordinal number of the electrode (30) which produced 
the sample numbered along said second direction, and wherein n is a function of a parameter k = i-g(j). 
where g is a function of i which relates the positions j of said electrodes in said array (1 2), numbered along 
said second direction, to the relative positions of said electrodes measured along said first direction in 

30 terms of L 



Patentanspruche 

35 1. Ein Verfahren fur das Erzeugen eines zweidimensionalen Bildes eines Charakteristikums einer Flache 
einer Wandung (14) einer Erdformation (16) unter Verwendung einer Matrix (12) von Sensoren (30), um- 
fassend eine Mehrzahl von Satzen von einem Oder mehreren Sensoren (30), wobei einige der Satze an 
unterschiediichen Positionen in der Matrix (12) angeordnet sind, gemessen lings einer ersten Richtung, 
und zu mind est einer der Sensoren (30) von mindestens einem der Satze unterschiedlich von alien in min- 

40 destens einem der anderen der Satze positioniert ist, gemessen langs einer zweiten Richtung senkrecht 

zu der ersten Richtung, wobei das Verfahren das Bewegen der Matrix (12) langs der Wandung (14) in 
einer ersten Richtung umfa&tzum Erzeugen eines Signal pegels an jedem Sensor (30), der reprasentativ 
ist fur den Wert des Charakteristikums in dem Abschnitt der Wandung (14) nahe zu ihm, Abtasten (36) 
derSignalpegel an alien Sensoren (30) der Satze bei jedem einer Aufei nan derfolge von Abtastzeitpunkten 

45 zum Erzeugen aufeinanderfolgender Originalrahmen von Sensorsignalabtastmustern, Verarbeiten (60) 

der Abtastmuster und Anlegen der Abtastmuster an ein zweidimensionales Wiedergabegerat (68) in Sat- 
zen derart, dad die verarbeiteten Abtastmuster von jedem der Originalrahmen Punkte in der Wiedergabe 
steuern, die in einem geometrischen Muster entsprechend dem geometrischen Muster der Satze von 
Sensoren (30) in der Matrix (12) bef indlich sind, gekennzeichnet durch, vor dem Anlegen der verarbeiteten 

so Abtastmuster an die Bildwiedergabe (68), Erzeugen (69) aus den Abtastmustern eines Korrektursignals 

n und Subtrahieren desselben von jedem der Abtastmuster zum Reduzieren irgendeines Artefakts in dem 
Bild, hervorgerufen durch Rauschen, das die Sensoren beeinfluBt, welches Rauschen jeden der Senso- 
ren (30) im wesentlichen gleich beeinfluBt wan rend jedes der Abtastzeitpunkte und in dem Bild einen Ar- 
tefakt erzeugt mit einer Form entsprechend dem Muster der Matrix (12) von Sensoren (30), wobei die Wer- 

55 te von n so erzeugt werden, dad abrupte Diskontinuitaten zwischen Werten jener Abtastmuster minimiert 

werden, die in unterschiediichen Originalrahmen von Sensoren (30) an derselben Position langs der er- 
sten Richtung erzeugt wurden, und wobei die Werte von n, die so subtrahiert werden, fur alle an demsel- 
ben Abtastzeitpunkt erzeugten Abtastmuster gleich sind. 
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2. Ein Verfahren nach Anspruch 1, bei dem das Korrektursignal n Werte aufweist, die ein Bilddiskontinui- 
tatsmaR I minimieren, worin I eine Funktion von Unterschieden in Signalabtast muster n ist, die uber eine 
Mehrzahl der Rahmen von Abtastmustern bewertet werden. 

3. Ein Verfahren nach Anspruch 1 oder2, beidemdieSensoren (J1-J27) derSatze, die an unterschiedlichen 
Positionen langs der ersten Richtung positioniert sind, langs dieser ersten Richtung Abstande (a) auf- 
weisen, die ganzzahlige Vielfache, einschlie&lich Eins, des Abstandes (d) sind, um welchen sich die Matrix 
(12) zwischen aufeinanderfolgenden Abtastintervallen bewegt 

4. Ein Verfahren nach Anspruch 2, bei dem I eine Funktion von S(i,j) - S(i,j+1) ist, worin S(i,j) = M(ij) - N(ij) 
ist, worin: 

i ein ganzzahliger Index ist, der die Position (a) jedes der Sensoren (30) reprasentiert relativ zu 
einer Referenz posit ion, wenn sein Signalpegel abgetastet wird, gemessen langs der ersten Richtung und 
in Ausdrucken von Einheiten gleich der Distanz (d), welchedie Matrix zwischen aufeinanderfolgenden Ab- 
tastungen durchlauft; 

j ein ganzzahliger Index ist, welcher die Position jedes der Sensoren (30) langs einer zweiten Rich- 
tung senkrecht zur ersten Richtung reprasentiert und gemessen ist in Ausdrucken der Ordnungszahl jeder 
Sensorposition, benummert langs der zweiten Richtung; 

M(ij) der gemessene Wert von jedem der Sensorsignalabtast werte ist, genommen an Sensorpo- 
sitionen (i j,); 

N(i,j) der Wert des Rauschens in dem Signalabtastmuster an jeder entsprechenden Sensorposition 
(ij) ist und im wesentlichen gleich ist fur alle Sensoren (30) zu irgendeinem gegebenen Abtastzeitpunkt; 
und 

S(i,j) der wahre Wert des Sensorsignalabtast musters ist, erzeugt von jedem der Sensoren an Po- 
sition (i,j.) bei Abwesenheit von Rauschen N(iJ). 

5. Ein Verfahren nach Anspruch 4, bei dem die Bilddiskontinuitatsfunktion I eine Funktion der absoluten 
Werte der Differenzen S(i,j) - S(i ist, bewertet uber eine Mehrzahl der Rahmen von Abtastmustern. 

6. Ein Verfahren nach Anspruch 4 Oder 5, bei dem die Bilddiskontinuitatsfunktion I eine Funktion derSumme 
der Quadrate der Differenzen S(ij) - S(i,j+1) ist, bewertet uber eine Mehrzahl der Rahmen von Abtast- 
mustern. 

7. Ein Verfahren nach Anspruchen 4 bis 6, bei dem das Korrektursignal n Werte aufweist, welche den Satz 
simultaner, linearer Gleichungen erfullen, ableitbar von djd n - 0. 

8. Ein Verfahren nach Anspruch 7, bei dem das Korrektursignal Werte besitzt gleich jenen, die aus einer 
Bewertung von: 

D 

n(k)= 2 T l (fcx).y (x) 
x-1 

resultieren, worin: 

n(k) der abgeschatzte Wert des Rauschens in jedem der Sensorsignalabtastmuster ist, gemessen 
wahrend desselben Rahmens k der Abtastung, 

x ein ganzzahliger Summationsparameter ist, der Werte von 1 bis D annimmt, 

D die Anzahl von Rahmen ist, welche im wesentlichen das gesamte Bild umfassen, 

T _1 (k, x) das (k, x)-Element der Matrix T~ 1 ist, 

T" 1 das Inverse einer gebandeten, symmetrischen Toeplitz-Matrix T reprasentiert, 
T die Matrixform der Koeff izienten der Rauschen enthaltenden Terme in dem Satz von linearen, 
simultanen Gleichungen reprasentiert, erhaltlich durch Setzen von d\fd n = 0, 

E-l 

y(x) = 2 (M(x+gOX j) • M(x+g(j). j+1) . M(x+g(j+l)J) + M(x+g(j+l), j+1)} 
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E die Anzahl der Sensoren (30) in der Sensormatrix (12) reprasentiert, und 
g(j) in Einheiten von i die Position von jedem der Sensoren (30) langs der ersten Richtung repra- 
sentiert relativ zu einer Referenzposition, die festliegt relativ zu der Matrix (30), ausgedruckt als eine 
Funktion der Position jedes Sensors (30) langs der zweiten Richtung (j). 

9. Ein Verfahren nach Anspruch 7, bei dem das Korrektursignal Werte aufweistgleich jenen, herruhrend von 
einer Bewertung von: 

w w 

n(k)=l/w I I T-^ki.k^^ydc+ki -k2) 
ki«l k2=I 



worin: 

n(k) der abgeschatzte Wert des Rauschens in jedem der Sensorsignalabtastmuster ist, gemessen 
wahrend desselben Rahmens k der Abtastung, 

w die GroSe des Filterfensters ist, die aqurvalent ist der Abmessung langs der ersten Richtung der 
GroGe des Abschnitts des Bildes, worin w irgendeine Zahl kleiner ats die Anzahl der Rahmen ist, welche 
die Gesamtheit des Bildes umfassen, 

ki und k 2 ganzzahlige Summations para meter sind, welche Werte von 1 bis w annehmen, 

"M^ka) das (k 1( k 2 )- Element der Matrix T~ 1 reprasentiert, 

T _1 das Inverse einer gebandeten, symmetrischen Toe pi itz- Matrix T reprasentiert, 

T in Matrixform die Koeff izienten der Rauschen enthaltenden Terme in dem Satz von linearen, si- 
multanen Gleichungen reprasentiert, die man erhSIt durch Setzen von djd n = 0, 

y(k) die Funktion der gemessenen Daten reprasentiert, wodurch: 



>(x) = Z (M(x+g(j)J) • M(x+g(j), j+1) - M(x+g(i+D, j) + M(x+g(j+l), 

30 j^l 



wobei E die Anzahl der Sensoren in der Sensormatrix reprasentiert, und 

g(j) in Einheiten von i die Position jedes der Sensoren (30) langs der ersten Richtung reprasentiert 

relativ zu einer Referenzposition, die festliegt relativ zu der Matrix (30), ausgedruckt als eine Funktion 

der Position jedes Sensors (30) langs der zweiten Richtung (j). 

10. Ein Verfahren nach Anspruch 9, bei dem die ganzzahlige Zahl 4 ist und die Funktion g(j) ist i-4R[j], worin 
RO] der Rest ist, resul tie rend von der Division von j durch 4. 

11. Ein Verfahren nach einem der vorangehenden Anspruche, bei dem die Sig nal a btast muster tiefenjustiert 
werden vor der Subtraktion des Korrektursignals n durch Auswahlen, aus den Rahmen von Abtastmu- 
stern, erzeugt durch Sensoren an unterschiedlichen Tiefenpositionen, korrespondierenderSatze von Ab- 
tastmustern, erzeugt von den Sensoren, wenn sie an derselben Tiefe sind. 

12. Ein System fur das Erzeugen eines zweidimensionalen Bildes mindestens eines Segments einer Ober- 
f I ache einer Erdformationswandung (14), umfassend: 

(a) eine zweidimensionale Matrix (12) von Sensorelektroden (30) und Mittel (27,26) fur das Bewegen 
der Matrix (12) langs des Wandungssegments (14) in einer ersten Richtung, wobei das Zentrum der 
Elektroden (30) von einander sowohl in der genannten ersten Richtung als auch in einer zweiten Rich- 
tung rechtwinkJig zu der ersten Richtung beabstandet ist, wobei mindestens einige der Elektroden (30) 
an unterschiedlichen Position en langs der ersten Richtung nicht ausgefluchtet miteinander langs der 
ersten Richtung sind; 

(b) Abtastmittel (36,56) f ur wiederholtes Abtasten der Signalpegel an alien Elektroden (30) zur Bildung 
entsprechender Originalrahmen von Signaiabtastmustern; 

(c) Tiefeneinstellmittel (60) fur die Verarbeitung der Originalrahmen von Signaiabtastmustern zum Er- 
zeugen von tiefenjustierten Rahmen der Signalabtastmuster, wobei alle Abtastmuster in jedem der tie- 
fenjustierten Rahmen an derselben Position langs der Wandung (14) produziert werden; 



30 



EP0 426 564 B1 



(d) eine zweidimensionale Wledergabeeinrichtung (60), die auf die tiefenjustierten Rahmen von 
Signalabtastmustern reagiertzum Erzeugen aufeinanderfolgender Zeilen in dem Bild in Reaktion auf 
entsprechende unterschiedliche der tiefenjustierten Rahmen von Signalabtastmustern, wodurch ein 
zweidimensionales Bild des Wandungssegments (14) gebildet wird, 
s gekennzeichnet durch: 

(e) Filtermittel (62) fur das Filtern der tiefenjustierten Signalabtast muster zum Reduzieren von Ar- 
tefakten, eingefuhrt in das Bild durch Rauschen, das alle Elektrodensignalpegel im wesentlichen im glei- 
chen MaBe beeinf lu&t, welche Artefakte den Mustern von Elektroden (30) in der Matrix (12) entsprechen, 
welche Filtermittel (62) Rauschsignalabschatzmittel umfassen, die auf die tiefenjustierten Abtastmuster 
10 reagieren zum Erzeugen von Korrektursignalen mit diskreten Werten n(k), identif iziert mit jedem der tie- 

fenjustierten Abtastmuster, und Mitt el fur das Subtrahieren jedes Wertes n(k) der Korrektursignale von 
dem entsprechenden Abtastmuster zum Reduzieren der Artefakte in dem Bild, und wobei die Korrektur- 
signalwerte ausgewahlt werden zum Minimieren von Differenzen zwischen den Werten der tiefenjustier- 
ten Signalabtastmuster, produziert durch langs der zweiten Richtung benachbarte Elektroden, und zum 
15 Verteilen der Korrektursignale gleichermaCen fiber Satze von tiefenjustierten Abtastmustern, erzeugt 

durch die gleichen Originalrahmen der Signalabtastmuster. 

13. Ein System nach Anspruch 12, bei dem die Werte N(k) so erzeugt werden, dap die Differenzen S(i,j) - 
S(i,j+1) minimiert werden, worin S die wahren Werte der tiefenjustierten Signalabtastmuster beim Fehlen 

20 von Rauschen darstellen, i ein Index ist, der die Tiefen reprasentiert, bei welchen die tiefenjustierten 

Signalabtastmuster ursprunglich erzeugt wurden, und j ein Index ist, der die Ordnungszahl der Elektrode 
(30) reprasentiert, welche das Abtastmuster erzeugte, benummert langs der zweiten Richtung, und wobei 
n eine Funktion eines Parameters k = i-g(j) ist, worin g eine Funktion von i ist, die die Posittonen j der 
Elektroden in der Matrix (12), benummert langs der zweiten Richtung, mit den relativen Position en der 

25 Elektroden, gemessen langs der ersten Richtung, in Ausdrucken von i in Beziehung setzt. 



Revendications 

30 1. Procede pour produire une image bidimensionnelle d'une caracteristique d'une zone de paroi (14) d'une 
couche de terrain (16) utilisant un ensemble (12) de capteurs (30) com pre n ant une plurality de jeux d'un 
ou de plusieurs capteurs (30), certains desd its jeux etant en des positions differentes dans ledit ensemble 
(12), tel que mesure le long de ladite premiere direction, et au moins Tun des capteurs (30) d'au moins 
Tun desdits jeux Etant positionne diffEremment par rapport a tous les capteurs dans au moins Tun desdits 

35 autres jeux, tel que mesure le long d'une second e direction transversal e a ladite premiere direction, ledit 

procede comprenant le deplacement de i'ensemble (12) le long de ladite paroi (14) dans une premiere 
direction, pour produire un niveau de signal a chaque capteur (30) representatif de la valeur de ladite ca- 
racteristique dans la partie de ladite paroi (14) adjacente correspondante, I'echantillonnage (36) desdits 
niveaux de signaux a tous lesdits capteurs (30) desdits jeux a chacune des successions de temps 

40 d'echantillonnage, pour produire des trames originales successives d'Echantillons de signaux des cap- 

teurs, le traitement (60) desdits Echantillons et Tap plication desdits Echantillons a un dispositif de visua- 
lisation bidimensionnel (68) dans des jeux, de maniere que lesdits Echantillons traites de chacun desdits 
points de contrdle de trames originales dud it dispositif de visualisation dont la configuration geometrique 
correspond a la configuration geometrique desdits jeux de capteurs (30) dans ledit ensemble (12), ca- 

45 racterise, avant 1'application desdits echantillons traites au dispositif de visualisation d'images (68), par 

la generation (69) a partir desdits Echantillons d'un signal de correction n, et la soustraction de celui-ci 
de chacun desdits Echantillons, pour red u ire tout artefact dans ladite image produite par le bruit affectant 
lesdits capteurs, ledit bruit affectant chacun desdits capteurs (30) de maniere sensiblementidentique pen- 
dant Tun quelconque des temps d'echantillonnage, et la production d'un artefact dans ladite image ayant 

so une forme correspondant a ladite configuration dudit ensemble (12) de capteurs (30), dans lesquels les 
valeurs de n sont generEes de maniere a reduire les differences brutales entre les valeurs de ces Echan- 
tillons produits dans differentes trames originales par des capteurs (30) a la meme position le long de 
ladite premiere direction, et dans lesquels les valeurs de n ainsi soustraites sont Egales pour tous les 
Echantillons produits au mEme temps d'echantillonnage. 

55 

2. ProcEdE selon la revendication 1 , dans lequel ledit signal de correction n a des valeurs qui reduisent une 
mesure de discontinue d'image I, ou I est une fonction des diffErences des Echantillons de signaux Eva- 
lues sur une plural itE desdites trames d'Echantillons. 
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Precede selon la revendication 1 ou 2, dans lequel lesdits capteurs (J1 a J27) desdits ensembles qui sont 
places en differentes positions te long de ladite premiere direction sont espaces le long de ladite premiere 
direction par des distance (a) qui sont des multiples, un compris, de la distance (d) dont se deplace ledit 
ensemble (12) entre des intervalles d'echantillonnage successifs. 

Proc6de selon la revendication 2, dans lequel I est une fonction de S(i,j)-S(i,J+1), ou S(iJ)= M(ij)-N(iJ), 
ou: 

i est un indice entier representant la position (a) de chacun desdits capteurs (30), par rapport a 
une position de reference, lorsque son niveau de signal est echantillonne, mesure le long de ladite pre- 
miere direction et en fonction d'unites egales a la distance (d) de deplacement dudit ensemble entre les 
echantillonnages successifs; 

j est un indice entier representant la position de chacun desdits capteurs (30) le long d'une seconde 
direction transversal e a ladite premiere direction, et mesure en fonction du nombre ordinal de chaque po- 
sition de capteurs, numerates le long de ladite seconde direction; 

M(i,j) est la valeur mesuree de chacun desdits echantillons de signaux de capteurs relevee aux 
positions des capteurs (i j); 

N(iJ) est la valeur dudit bruit dans I'echantillon de signal a chaque position de capteurs correspon- 
dante (ij), qui est essentiellement la meme pour tous lesdits capteurs (30) a un temps d'echantillonnage 
quelconque donne; et 

S(i,j) est la valeur vraie de I'echantillon de signal des capteurs produit par chacun desdits capteurs 
a la position (ij), en I'absence dudit bruit N(i,j). 

Proc6de selon la revendication 4, dans lequel ladite fonction de discontinuity d' image I est une fonction 
des valeurs absolues desdites differences S(i,j)-S(i,j+1), evalu£e sur une plural ite desdites trames 
d'echantillons. 

Proc6d6 selon la revendication 4 ou 5, dans lequel ladite fonction de discontinuity d' image I est une fonc- 
tion de la somme des carres desdites differences S(i,j)-S(i,j+1), evalu£e sur une plurality desdites trames 
d'echantillons. 

Proc^de selon les revendi cat ions 4 a 6, dans lequel ledit signal de correction n a des valeurs qui satisfont 
I'ensemble des equations lineaires simultanees derivables de d\/d n =0. 

Precede selon la revendication 7, dans lequel ledit signal de correction a des valeurs egales a cedes re- 
sultant d'une evaluation de: 

D 

n(k) = I t-1 (k, x) • y (x), 

x=l 

ou: 

n(k) est la valeur estimee dudit bruit dans chacun desdits echantillons de signaux de capteurs me- 
suree pendant ladite trame k identique d'echantillonnage, 

x est un parametre de sommation entier qui peut prendre des valeurs comprises entre I et D, 

D est le nombre de trames comprenant essentiellement I'image entiere, 

T- 1 (k. x) represente I'6l6ment (k, x) de la matrice T- 1 , 

T-i represente Tinverse d'une matrice bande symetrique de Toeplitz T, 

T represente sous forme de matrice les coefficients des termes contenant le bruit dans I'ensemble 
des equations simultanees lineaires susceptibles d'dtre obtenu en posant d\fd n = 0. 

E-l 

y(x)=I{M<x+g(j) r j>-M(x+g(j) f j + l)- M (x+g<j + l),j )+ M(x+g(j + l,,j + i )> 

E represente le nombre desdits capteurs (30) dans ledit ensemble de capteurs (12), et 
g(j) represente en unites de i la position de chacun desdits capteurs (30) le long de ladite premiere 
direction par rapport a une position de reference determinee en relation audit ensemble (30), exprimee 
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en tant que fonction de la position de chaque capteur (30) le long de ladite seconds direction (j)- 

9. Procede selon la revendication 7, dans lequel ledit signal de correction a des valeurs egales a ceiles re- 
sultant d'une evaluation de: 

5 

w w 

n(k) = 1/w I I T-l{k 1# k 2 ) • y (Jc+k! - k 2 ) 
k 2 =l k 2 »l 

10 

ou: 

n(k) est la valeur estimee dudit bruit dans chacun desdits echantillons de signaux de capteurs me- 
sur6e pendant ladite trame k identique d'echantillonnage, 

w represente la taille de la fenetre de f iltre qui est equivalente a la dimension le long de ladite pre- 
15 miere direction de la taille de ladite partie de ladite image, ou w est un nombre quelconque inferieur au 

nombre de trames que comprend la total ite de ('image, 

et k 2 represented les parametres de sommation entiers qui peuvent prendre des valeurs compri- 
ses entre 1 et w, 

T^ki.kJ) represente le (k^kj element de la matrice T-\ 
20 I -1 represente Tin verse d'une matrice bande symetrique de Toeplitz T, 

T represente sous forme de matrice les coefficients des termes contenant le bruit dans ('ensemble 
des equations simultanees lineaires susceptibles d'etre obtenu en posant d\id n = 0. 
y(k) represente la fonction desdites donnees mesurees, grace a quoi: 



25 

E-i 

y(x) S L { M(x + g(j),j ) .H(x + g(j ) ,j + l)-M(x + g(j + l), j >+M<x+g< j+ l > , j+ i> > 
j = l 

30 

E represente le nombre desdits capteurs dans ledit ensemble de capteurs, et 
g(j) represente en unites de i la position de chacun desdits capteurs (30) le long de ladite premiere 
direction par rapport a une position de reference determinee en relation audit ensemble (30), exprimee 
en tant que fonction de la position de chaque capteur (30) le long de ladite seconde direction (j). 

35 

10. Procede selon la revendication 9, dans lequel ledit nombre entier est 4 et ladite fonction g(j) est i-4R[j], 
ou R[j] est le reste resultant de la division de j par 4. 



11. Procede selon I'une des revendications precedentes, dans lequel lesdits echantillons de signaux sont 
ajustes en profondeur avant soustraction dudit signal de correction n en selection nant, a part ir des trames 
d'echantillons produits par des capteurs a differentes profondeurs, des jeux correspondents d'echantil- 
lons produits par lesdits capteurs lorsqu'ils se trouvent a des profondeurs identiques. 



12. Systeme pour produire une image bidimensionnelle d'au moins un segment d'une surface d'une paroi de 
couche de terrain (14), comprenant: 

(a) un ensemble bidimensionnel (12) d'electrodes de capteurs (30) et des moyens (27, 26) pour de- 
placer ledit ensemble (12) le long dudit segment de paroi (14) dans une premiere direction, le centre 
desdites electrodes (30) etant espace I'une de I'autre, a la fois, le long de ladite premiere direction et 
le long d'une seconde direction a angles droits par rapport a ladite premiere direction, au moins cer- 
taines desdites electrodes (30) a differentes positions le long de ladite premiere direction etant non 
alignees I'une par rapport a I'autre le long de ladite premiere direction; 

(b) des moyens d'echantillonnage (36, (56) pour echantillonner de maniere repetitive les niveaux de 
signaux a touts s lesdites electrodes (30) pour former des trames originates correspond antes d'echan- 
tillons de signaux; 

(c) des moyens d'ajustement en profondeur (60) pour t ratter lesdites trames originates d'echantillons 
de signaux pour produire des trames ajustees en profondeur desdits echantillons de signaux I'ensem- 
ble desdits echantillons de chacune desdites trames ajustees en profondeur etant prod u it a la meme 
position le long de ladite paroi (14); 
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(d) un dispositif de visualisation bidimensionnel (60) sensible auxdites trames ajustees en profondeur 
d'echantiiions de signaux pour produire des lignes success ives dans ladite image en reponse aux tra- 
mes differentes correspondantes desdites trames ajustees en profondeur d'echantiiions de signaux, 
pour former ainsi une image bidimensionnelle dudit segment de paroi (14), 
caracterise pan 

(e) un moyen de f iltrage (62) pour f iltrer lesdits echantillons de signaux ajustes en profondeur pour 
reduire les artefacts introduits dans ladite image par le bruit qui affecte tous lesdits niveaux de signaux 
d'electrodes essentiellement de maniere identique, lesdits artefacts correspondant aux configurations 
d'electrodes (30) de I'ensemble (12), ledit moyen de f iltrage (62) comprenant un moyen d'estimation du 
signal de bruit sensible auxdits echantillons ajustes en profondeur pour generer des signaux de correction 
ayant des valeurs discretes n(k) identifies avec chacun des echantillons ajustes en profondeur, et un 
moyen pour soustraire chaque valeur n(k) desdits signaux de correction de I'echantillon correspondant 
pour reduire lesdits artefacts dans ladite image, et dans lequel lesdites valeurs de signaux de correction 
sont selectionnees pour reduire les differences entre les valeurs desdits echantillons de signaux ajustes 
en profondeur produits par des electrodes adjacentes les unes aux autres le long de ladite seconds di- 
rection, et pour re part ir lesdits signaux de correction de maniere egale sur les jeux d'echantiiions ajustes 
en profondeur produits par les memes trames originates desdits echantillons de signaux. 

Systeme selon la revendication 12, dans lequel lesdites valeurs N(k) sont generees pour reduire les dif- 
ferences S(i,j)-S(ij+1 ), ou S represente les valeurs vraies desdits echantillons de signaux ajustes en pro- 
fondeur en I'absence dudit bruit, i est un indice representant les profondeurs auxquelles lesdits echantil- 
lons de signaux ajustes en profondeur ont ete produits originalement, et j est un indice representant le 
nombre ordinal de ('electrode (30) qui a produit I'echantillon numerate le long de ladite seconde direction, 
et dans lequel n est une fonction d'un para met re k = i-g(j), ou g est une fonction de i qui etablit un rapport 
entre les positions j desdites electrodes dudit ensemble (12), numerotees le long de ladite seconde di- 
rection, et les positions relatives desdites electrodes mesurees le long de ladite premiere direction en 
fonction de i. 
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SAMPLE VALUES BEFORE DEPTH ADJUSTMENT 



.5 
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3 
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7 




7 


7 


\ 7 


47 


- 7 


7 
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9 


9 


9 


9 


9 
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9 
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10 


10 


10 
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11 
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12 


12 
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13 
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13 
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14 
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14 
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DEFINE g ( j ) 



134 



SUBSTITUTE IN 



DEFINE E 



136 
SUBSTITUTE 



SUBSTITUTE IN 



E-1 



V 1 |2n(k)-n[k+g(j)-gG+1)] 
J=1 

- n[k+g(j+1)-g( j )]} = y( k ) 



IN 



130 



A 



32 



138 



E-1 



y( K ) = E |M(k+g( j ) , j ) - M(k+g( j ), j+1) - 
j =1 M(k+g0+1) . j ) + M(k+g0+1), j+1)} 



160 



164. 



EXPAND FOR 
j=1 TO (E-1) 
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EVALUATE FOR 
g( j ) AND gG+1) 



EXPAND FOR 
j=1 TO (E-1) 



ACCESS 

M( i . j ) 

FOR 
i=1 TO D 
j=1 TO E 
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COMBINE TERMS AND SIMPUFY 
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GENERATE 
Hy( k ) FOR 
j=1 TO D 



r 



DEFINE D 



162 



■146 



144— 



GENERATE MATRIX T 
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DEFINE W 



170 SUBSTITUTE 



IN 



— DEFINE (» x w) 
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W W 

n(k) = 1/w J] J jT~ (k k) y(k+k -K )} 
k 1= 1 k 2 =1 12 12" 



152 



2l 



1 



GENERATE 
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-1 



EXPAND FOR 
k!= 1 TO w 

k 2 = 1 TO w 
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GENERATE n( k ) 
FOR k=1 TO D 

1 



178 



S( i , j ) = M( i , j ) - n( k ) 
FOR i-1 TO D; j=1 TO E 



180 



SAVE S( i, j ) FOR IMAGE 



FIG. 10 
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FIG . 11 
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